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We present a theory of the elasticity and fluctuations of the Smectic A and C phases in uni- 
axial, anisotropic disordered environments, e.g., stretched aerogel. We find that, bizarrely, the 
low-temperature, lower-symmetry Smectic C phase is less translationally ordered than the high- 
temperature, higher-symmetry Smectic A phase, with short-ranged "m = 1 Bragg glass" and alge- 
braic "XY Bragg glass" order, respectively. The AC phase transition belongs to a new universality 
class, whose fixed points and exponents we find in a d = 5 — e expansion. We give very detailed 
predictions for the very rich light scattering behavior of both phases, and the critical point. 
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I. INTRODUCTION 

Randomly pinned elastic media occur in many con- 
texts, including disordered superconductors Jl| , charge 
density waves Josephson junction arrays |3j and He- 
lium in aerogel p\. By far the most exotic phenomenon 
(at least in the - admittedly biased - opinion of the cur- 
rent authors) in such systems occurs in liquid crystals 
in aerogel: anomalous elasticity. That is, many liquid 
crystals in aerogel exhibit scalings of their elastic ener- 
gies that differ radically (specifically, by non-trivial power 
laws) from those found in the absence of pinning. 

However, despite the considerable amount of prior 
work done on these problems, there had, until recently 
0, been no previous theoretical work on phase transi- 
tions in pinned liquid crystal systems and little on phase 
transitions in any pinned elastic system. In this paper 
we remedy this by treating the smectic A to smectic C 
(hereafter AC ) transition in an anisotropic, uniaxial dis- 
ordered environment. A brief summary of a few of our 
results has already appeared in Q . Such an environment 
could be realized, e.g., by absorbing the liquid crystal in 
uniaxially stretched aerogel. For brevity, we will here- 
after refer to the special uniaxial direction as the z-axis 
of our co-ordinate system. 

The AC phase transition in such a system separates the 
two novel, anomalously elastic glassy phases discovered 
(in a totally different context) and treated in reference[8|. 
The high temperature phase (T > Tac) is the glassy 
analog of the smectic A phase of the pure problem, in 
that both the layer normal N and nematic director n lie, 
on average, along the z-axis. This phase is in the uni- 
versality class of the random field XY modeli4]; hence, 
following P, we call it the "random field XY smectic 
Bragg glass", or XYBC for short. 

The low temperature phase is the glassy analog of the 
smectic C phase, in that both the layer normal N and 
nematic director n tilt from z-axis. This tilting obeys 
the "geometrical constraint" that N , n and z are in the 
same plane, with z between TV and h and the angles be- 
tween these three vectors non-fluctuating, albeit temper- 



ature dependent. Hence, as in the smectic C phase in an 
isotropic environment, the only new "Goldstone mode" 
degree of freedom associated with the tilting is the over- 
all azimuthal angle of rotation of the vectors N and n 
about z. 

The temperature dependence of the angles 9l{T) be- 
tween N and z and 9n{T) between fi and z are the same 
near T^, that is, the ratio BlI^ti is a (negative) constant 
near T^- Hence, we are free to choose either the layer 
normal tilting angle (T) , or the molecular tilting angle 
9n(T), both of which can be measured experimentally, as 
the magnitude of the order parameter for the AC transi- 
tion. 

The smectic C phase is in the universality class of the 
"to = 1 smectic Bragg glass" phase studied in [8j] . 

We call both the A and C phases "glassy" because 
both lack long-ranged translational order due to the dis- 
ordering effect of the random environment (i.e. the aero- 
gel). The extent of this destruction, however, differs 
greatly between the two phases. Strikingly, it is the low- 
temperature, higher-symmetry, "Smectic C glass" that 
has less translational order. In the "glassy smectic A' or 
"XYBC " phase, translational correlations are "quasi- 
long-ranged" , by which we mean they decay as power 
laws with distance. In the "glassy smectic C" or "to = 1 
Bragg glass" phase, these correlations are short-ranged. 
These differences in the translational correlations lead to 
radically different X-ray scattering signatures in the two 
phases which we will now describe. 

In the "glassy A" or "XYBC" phase, the X-ray scat- 
tering intensity I{q} diverges near the smectic Bragg 
peaks, which occur at g = nqoz for all n is integer, where 
qa = 27r/a, with a the smectic layer spacing. Because 
of the lack of true, long-ranged translational order, this 
divergence is not in the form of a delta- function; rather, 
it is an integrable power-law divergence: 



aql] 



(1) 



where a is a non-universal constant of order 1 and q± is 
the magnitude of the projection of cj perpendicular to z. 
Note that the power law —3-1- .55n^ characterizing the 
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divergence of the nth peak depends on which peak we are 
considering. Indeed, only the first 2 peaks {n ~ 1 and 
n = 2) actually diverge. 

In contrast, in the "glassy C" or "m = 1 _BG" phase, 
the peaks in the X-ray scattering intensity are broad, 
with I{q) finite for all q. 

As T ^ Tag from above (i.e., on the A side), the 
sharpness of the peak disappears in an unusual way. 
The peaks look broad and qualitatively Lorentzian for 
q 's sufficiently far from the Bragg peak position nqoz, 
while for q 's sufficiently close to the peak, it follows 
the power law divergence (H]). "Sufficiently close", in 
this context, means that both \q±\ ^ Sq']_{n,T), and 
\qz — nqo\ ^ Sq^iT), where the crossover wavevectors 
Sqlin.T), Sqlin.T) are given by 

Sqlin.T) « , (2) 

5ql{n,T) (X {iiy^^^ , (3) 



l(q,,q =fO) 




FIG. 1: The gz-dependence of the X-ray scattering intensity 
for gx = in the smectic A phase. In the C phase, the sharp, 
power law peaks disappear, leaving only the broad peak. 



where n again is the index of the Bragg peak, and and 
are correlation lengths along and perpendicular to the 
smectic layers respectively that both diverge extremely 
strongly as T ^ Tag- Specifically: 



^1., cx exp(^ |T - 



TagV""), 



(4) 



where is a universal exponent which we have calcu- 
lated in the e-expansion discussed below, and A is a non- 
universal constant. 

These X-ray scattering predictions are illustrated in 
Fig. [1] and [2l The divergence of z implies that, as 
T T^fj, the algebraic "spikes" on the top of the broad 
quasi-long-ranged peaks get narrower and less intense, 
vanishing completely at Tag- Lowering temperature fur- 
ther leads only to the broad peaks of the Smectic C phase. 
This entire scenario of sharp peaks at high temperature 
and broad peaks at low temperature is completely coun- 
terintuitive, and contrary to the behavior seen in almost 
every other translationally ordered system [l^ . 

We turn now to the smectic C phase, and in particular 
to the question of why this phase, though of lower sym- 
metry than the smectic A phase, is less translationally 
ordered. In fact, it is precisely the new broken symmetry 
of the smectic C phase-that is, the tilt of the layer normal 
and nematic director-that causes this. This is because, 
while the energetically preferred layer normal N and ne- 
matic director n in the smectic A phase are unique-they 
must point along z-there are infinitely many energetically 
preferred orientations of N and h in the C phase: N can 
lie anywhere on a cone making an angle Ol{T) with z, 
which, combined with the "geometrical constraint" (de- 
scribed in the third paragraph) on the directions of N 
and n, also determines the similar cone of directions h 
can point. This exact symmetry of the elastic energy of 
the smectic C phase means that the direction perpendic- 
ular to the z— < N > plane (where < iV > is the mean of 
N; i.e., the direction of spontaneous tilt) becomes "soft": 



l(q,=q„) 


^l(q)~|qf" 







FIG. 2: The gx-dependence of the X-ray scattering intensity 
for Qz = qo in the smectic A phase. Again, the sharp peak 
vanishes in the C phase. 



that is, an easy direction for layer displacements to vary 
in. Precisely such softness occurs (for different reasons) 
in the "m = 1 smectic" studied in Q and, indeed, ex- 
panding our elastic Hamiltonian around the new broken 
symmetry, tilted smectic C ground state, we obtain an 
elastic Hamiltonian describing small positional fluctua- 
tions about the tilted smectic C ground state, which, 
after some manipulation, proves to be identical to that 
studied for "to — I smectic" in Q. Thus, we can simply 
transcribe the results of 8] to this problem. 

The most striking result of reference [sj is that this 
phase exhibits a bizarre phenomenon known as "anoma- 
lous elasticity" . That is, the elastic moduli that char- 
acterize the elastic response of this phase can no longer 
be called "elastic constants" , because they are not con- 
stants. Rather, they become singular functions of the 
length-scale, or, equivalently, the wavevector, at which 
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they are measured. 

The elastic constants that characterize the sniectic C 
phase in uniaxial disordered media are defined in (|206p . 
which gives the "replicated" elastic Hamiltonian describ- 
ing the disordered, glassy smectic C phase. Of these, 
the bend modulus and the disorder variances or ef- 



where the wavevector- independent quantities Kc(T), 
A^,_^,(T), and 7c(T) are the "half-dressed" (i.e., renor- 
malized by the critical fluctuations we will discuss in a 
moment, but unrenormalized by smectic C fluctuations) 
values of the corresponding elastic moduli and disorder 
variances; we've defined a new, non- orthogonal set of 
wavevector coordinates: 

qx- = qx- r^z, (8) 
qs' = qs, (9) 
qz' = qz (10) 

where qs is the component of q perpendicular to the z- 
< N > plane, and qx is the component of q within the 
z-< N > plane but orthogonal to z, £,±{T) and ^z{T) are 
temperature-dependent lengths, and T(T) is a dimension- 
less temperature-dependent constant. The temperature 
dependences of ^^(r), ^.(T), K,{T), A-,^^,{T), 7,(T) 
and r(T) are all singular at the transition temperature 



fective disorder strengths Ag' and A^' diverge, while the 
in-plane tilt modulus 7 (which gives the energy cost for 
changing the angle between N and z from 9l{T)) van- 
ishes, as \q\ 0, according to 



(5) 



(6) 



(7) 



I 

Tac, with a(r), UT), K,[T), A^,^,,(r), and V{T) di- 
verging, and jdT) vanishing, as T ^ Tac- These scaling 
laws governing their respective divergences and vanish- 
ings will be given later in this introduction, in section 
IVIIIl and in Appendix O 



Here and throughout the rest of this paper, a tilde 
{~ ) over an exponent denotes an exponent describing 
anomalous elasticity of the C phase. 



The Mm?;ersaZ exponents fiK,f}-y,fis' ,x' ,Cx' and Cz' ap- 
pearing in the above expressions were calculated in [8|, 
using an approach that yielded numerical estimates and 
error bars for them. We will review the logic of their 
approach in Appendix \^ here, we merely quote their 



k,{T){Uqs'r^'' , [i^qs'P » Uqx',{i^qs'P » ^zqz- 
KcmiUqx')-^"/'^^', Uqx' » iC±qs')'^^',Uqx' » i^zqz')^ 
Kcim^zqz'r'""^^^', ^zqz' » iUqs'P.^zqz' » {Uqx')^ 



( A^,_,,(T)(^i(?,0"''='-', iUqs'P » Uqx', {^±qs'P » izqz' 

As',.' = I A^,^,,(r)(e±g:.')-'^='-'/^^^', ^^qx' » {Uqs'P,aqx' » (e.«.')_^ 

[ A^,_,,(r)(e.g.')-'^='-'/^^', ^zqz' » {Uqs')'^'',izqz' » {^±qx')^ 



f ic{T)iUqs')''-' , i^qs'p » Uqx',{i±qs'P » f.?.' 

lc{T){Uqx'r-'^^^', Uqx' » iUqs')^^' ,S.^qx' » {Lqz'P 
{ lc{T){Lqz'T-'~^'\ izqz' » i^±qs')^^',Lqz' » (Cife')^ 
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results: 



VK 


= 0.50 ±0.03, 


(11) 


V-t 


= 0.26 ±0.12, 


(12) 


Vs' 


= 0.132 ± 0.002, 


(13) 




= 1.62 ±0.08, 


(14) 


C.' 


= 1.75 ±0.02. 


(15) 



These exponents also obey the exact scaling relations 
(in d = 3): 

C.' = 2-(^), (16) 

= 2-^, (17) 

77,. = !h. + 2fjK-l- (18) 

In addition to the above results derived in reference Q, 
we derive in section IIIII an additional exact scaling rela- 
tion for fj.j.i: 

fj^, = 2 + fjs> - fiK - V-y , (19) 

from which it follows, upon inserting the numerical re- 
suhs Eqns. (HH [H [HI) that 

= 1.37 ±0.15 . (20) 

Having described the A and C phases in this uniaxially 
disordered system, we now turn to the transition between 
them. We find that the upper critical dimension due for 
this transition, below which it is no longer accurately de- 
scribed by a purely Gaussian theory, is due = 5. We have 
studied this phase transition in an e = 5 — d expansion, 
where d is the dimension of the space filled by the liquid 
crystal (i.e., d = 3 is the case of physical interest), and 
find that there is a stable renormalization-group fixed 
point, the existence of which implies a second-order phase 
transition with universal critical behavior. Of course, 
since due = 5, e = 2 in the physically interesting case 
d = 3, which is far too large an e for the results for ex- 
ponents that we quote to be quantitatively reliable in 
d = 3. However, we expect the qualitative features of 
the transition that we find to be robust down to d = 3 
(although we will discuss some rather strong caveats to 
this statement later). 

The second-order nature of the transition is, as usual, 
manifested in universal power-law dependence of many 
physical observables on reduced temperature T — Tac- 
In particular, the layer normal and director tilt angle 
OlAT) obeys 

0L,n{T) = Al.u {Tag - Tf , (21) 



where the order parameter exponent P is universal, and 
given, to leading order in e = 5 — d, by 

(^=\-T0+O{e'), (22) 
the amplitude Ai_„ are non- universal and depend on the 
system. 

We have also calculated the specific heat exponent 

^ = -^.+0{e')- (23) 



The order parameter Nj_ for this transition is the pro- 
jection of the smectic layer normals N onto the xy plane 
(i.e., perpendicular to the direction along which the aero- 
gel is stretched ). Above Tag, the two point real space 
correlations of N± decay rapidly with distance, with cor- 
relation lengths ^2 and parallel and perpendicular 
to z-axis respectively which behave quite differently as 
T 2^Ac" Specifically, both diverge as power laws in 
(T-Tag)- 

a.. « IT-TacI-""- (24) 

but with exponents v± and z/j, that differ from each other 
due to the strong anisotropy of the problem. From our 
e-expansion, we find 

For T < Tag, i-e., in the C phase, the correlation 
lengths and can still be defined, now by looking 
at the connected correlations of A^. They continue to be 
given by the scaling law (|24p . with the same exponents 
([25]) and 

In addition to their role as correlation lengths for N, 
the lengths and are also the ones that appear in 
the scaling laws Eqns. (O, (O and ([7]) for the anomalous 
elasticity in the C phase. 

We also find that, as in the C phase, the system ex- 
hibits anomalous elasticity right at Tag as well. Specif- 
ically, right at Tag, the smectic layer bend modulus 
K becomes strongly wavevector-dependent, vanishing as 
according to the scaling laws 
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IK 



-Vk 



(27) 



where rjK < 0, and the anisotropy exponent 



are given as 



This anisotropy exponent also obeys 



(28) 



(29) 



In Eq. ([77]) . the constant Kq is the "completely bare" 
(i.e., unrenormalized by any fluctuations) value of the 
bend modulus (i.e., the value it would have in the absence 
of the aerogel). While weakly temperature dependent, 
like all elastic moduli, Kq has no critical singularity near 
Tag- 

Also in Eq. (|77|) ^^'^ are the length scales in the _L and 
z directions beyond which the elasticity becomes anoma- 
lous. These lengths depend on the disorder strength, and 
hence the aerogel density, and remain finite (and non- 
zero) as T ^ Tac, although both diverge as the aerogel 
density pA 0. In three dimensions these length scales 



J 



647r Kl 



1/2 



5/4 



l/4^0l/2 



B, 



(30) 
(31) 



where is the " completely bare" value of disorder vari- 
ance A(, which is discussed below. 

The disordering effect of the random aerogel matrix 
can be quantified by disorder variances Aj and A^, de- 
scribing tilt and compressive stresses respectively. The 
former arise due to random torques exerted by the aerogel 
on the smectic layers, causing them to tilt, while the lat- 
ter are caused by random compression (or stretching) of 
the smectic layer spacing induced by the aerogel, causing 
layers to move closer together (or further apart). These 
variances also become anomalous at T = Taci obeying 



At_c (g,r = Tac 



AL(e^L9±) 



r 



-Vt,, 



(32) 



where A^ is the "complete bare" value of Ac. 

Because our model at T = Tac proves remarkably sim- 
ilar (although not quite identical) to that for a smectic A 
in isotropic disordered media, we expected, when we be- 



J 



B (g, T = Tac) = BoiCl^Lq^r^ ft 



gan our study, to find anomalous behavior for the smectic 
layer compression modulus B, with it vanishing as (f — s- 
according to 



(33) 



where Bq is the "complete bare" value of i?, as found [6| 
in the isotropic disordered smectic A problem. To our 
surprise, we found that, in this problem, B exhibits no 
such anomaly, remaining constant as g— > 0. This result 
is exacts not an artifact of the e = 5 — d expansion, and 
we expect it to hold in c? = 3. fK,t.c.B{x) are universal 



scaling functions, which have the property: 



fB,K,t,c{^) 



a; < 1 
, a; > 1 



(34) 



The exponents 'qK, % and r]c are non-zero, however. 
Both rjK and rjt are zero to 0(e), but non-zero to O(e^), 
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FIG. 3; Illustration of the three distinct regions in wavevector 
space with different wavevector-dependences of At^c and if, 
which is given in Eq. ([38}. 



and given by 

m = CKe' + 0(e3), (35) 
m = CAe' + 0(e3) (36) 

with Ck = (32 In (4/3) - 10) /225 -.00353 and Ca = 
(l21n(4/3) - i) /225 .01386. Note that Ck < 0, 
which implies that K vanishes as g — > 0. This is an- 
other unexpected result: K diverges as g — > in every 
other problem of this type [1, lisjpreviouslv studied. 
Of course, whether rjK remains negative all the way down 
from d = 5 to d = 3 remains an open question. 
The remaining exponent 

r?c-2-| + 0(a. (37) 
5 

For T bigger than Tag the (wavevector q and tem- 
perature T-dependent) disorder variance (q, T) and 
layer bend modulus K (q, T) are given by their T = Tag 
forms equations (|27p and ((32)) if either q±^± ^ 1 or 
9zCz ^ 1- Otherwise (i.e., if both q±£,± ^ 1 and 
« 1), At,c (X cc {T-TacT"^'''-' and K cx 
(x (T — Tc)~'^^''^', results which can be (and, in fact, 
were) straightforwardly obtained by matching the behav- 
ior for either q±^± ^ 1 or q^^z ^ 1 to the wavevector- 
independent behavior expected in the A phase when both 
q_i_(,±_ <^ 1 and q^^^ ^ 1- This can be summarized as fol- 
lows: 



At,,{q,T), i^(<f,r)« <^ q-"*-- region 2 , (38) 
[ e± region 3 

where (x {T — Tag) '^^ , and the three regions in q 
space are illustrated in Fig. [3] 

Having defined all of the critical exponents associated 
with this problem, we can now give several exact scaling 



relations between these exponents: 

a - 2-i.i(d-l + ^- rjt) , (39) 

P = iy^i2d-6 + 3rjK-2rjt)/4, (40) 

ly, = Ci^±- (41) 

Note that a does not obey the usual anisotropic hyper- 
scaling relation a = 2 — (d— l)i'± — Vz', this is due to the 
strongly relevant disorder. 

All of these exponents can be deduced experimentally 
from measurements cither for T > Tag or T < Tag- The 
specific heat can, of course, be measured by the usual 
thermodynamic measurements on either side of the tran- 
sition. 

The spontaneous tilt angle 9l(T) is another matter. 
In pure systems, this is simply the angle between the 
sharp Bragg spot of the smectic C and 5-axis. How- 
ever, as discussed earlier, in the Smectic C glass phase 
the X-ray scattering peak is broad, with a width that re- 
mains finite as T ^ ^^ac- Hence, close to Tag, the broad 
peak strongly overlaps z-axis, rendering measurement of 
Ol{T) by this approach impossible. Unfortunately this 
is precisely the temperature range one needs to study to 
measure (3. 

Fortunately, an alternative measure of Ol{T) is avail- 
able in the disorder averaged mean value of dielectric or 
diamagnetic susceptibility tensors Xij and e^. In the A 
phase, one of the principal axes of both tensors is along 
z-axis. In the C phase, this axis spontaneously rotates 
away from z-axis due to the spontaneous tipping of the 
smectic layers and nematic directors. The rotation angle 
is proportional to Ol{T). 

The best experimental probe of the critical phenomena 
(in particular, the correlation lengths and ^z), as well 
as of the anomalous elasticity and disorder-induced fluc- 
tuations in both the A and the C phases, is light scatter- 
ing, which probes fluctuations in both the dielectric(eij) 
and diamagnetic (xij) susceptibility tensors. 

Specifically,the large biaxial fluctuations on both sides 
of the transition lead to fluctuations in both and Xij 

that are proportional to [ll| (g) = {N^ {q) Nj- {-q}), 

where N± is the order parameter, which we remind the 
reader is the projection of the local layer normal perpen- 
dicular to the z-axis. 

The behavior of is radically different in the two 
phases, and quite involved (as the intrepid reader is about 
to discover) in both. Here we'll begin by summarizing the 
behavior in the (relatively) simpler A phase, and then 
proceed to a description of the even more complicated 
behavior in the C phase. 

In the A phase (T > Tag), there are three important 
contributions to Cij (q) : two arise from the previously 
discussed random tilt {At{q)) and compression (Ac(g)) 
disorders. The third is caused by random positional pin- 
ning Ap, which reflects the tendency of the aerogel to 
pin the smectic layers in particular positions. This type 
of disorder is irrelevant for T < Tag, but becomes im- 
portant at sufficiently long distance (or small q) in the 
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A phase (i.e., above Tag)- Indeed, as we will see in sec- 
tion IVIIIi this random positional pinning Ap actually 
dominates the fluctuations at small q in the A phase; 
furthermore, it is precisely this pinning that leads to the 
destruction of long-ranged translational order, and the 



replacement of sharp, delta-function Bragg peaks with 
power-law peaks, in the A phase. 

Combining the contributions from these three distinct 
types of disorder gives 



At (g, T) qi , A, (q, T) qlql 

,T)qif 

(42) 



Cijiq) — LjJq) Tj -I TT 

' V (Bq'z + D{T)ql + K {q, T) qlf {Bql + D{T)ql + K (g, T) qif 

CB^Dql \ 



3 



[Bql + D{T)ql + K (q, T) q^Y ql ) 

I 

where L^{q) = q:^qf /q\ is the projection operator along ature dependence of in each of these regions is: 
g^L J the compression modulus B remains constant as g* — *■ 
0, C is a universal, 0(1) constant (C « I.IOtt^), 



r {T~TAc)qz''°'^ region 1 
D{T) ^ I [T ^ TAc)qT" region 2 , (43) 
[ region 3 



where rjD = '2 — tjk — and the (wavevector q and 

temperature T-dependent) disorder variances At (q,T), 
Ac (g, T) and layer bend modulus K (g, T) are given by 
Eq. ([55]) . The C-term comes from the periodic random 
pinning potential. The three regions in (|43p is again il- 
lustrated in Fig. [3l 



Clearly, the expression (|42|) for is somewhat com- 
plicated, due to the fact that three distinct types of disor- 
der contribute to layer normal fluctuations. The situation 
is further exacerbated by the strong wavevector and tem- 
perature dependence of the elastic constants D and K. 
Fortunately, there are well-deflned regions of wavevector 
q, in each of which one of the three terms (i.e., the At, 
the Ac, or Ap-embodied in the C-term) dominate C^j, 
thereby simplifying the expression for it. Each of these 
regions can be further subdivided according to which of 
the terms Bq^, D{T)q\, and K{q,T)q\ dominates the 
common denominator Bql +D{T)q\+K (g, T) q\ of aU 
three terms. In addition, K{q,T), Ac,t{q,T) and D{T) 
themselves have the crossovers embodied in Eqs. (|38l 
which further subdivide some of the regions of g- 
space into subregions of distinct behavior. 

Painstakingly, but essentially straightforwardly, sort- 
ing out these different regions leads to the identification 
of eight (count 'em!) distinct regions of qualitatively dif- 
ferent wavevector-dependence for C^j, which are illus- 
trated in Fig. 3) The leading wavevector and temper- 



>'4i \ i± J 9J 



q± ' 

s 2-J7K 2 

± \2riK-'nt-i 



A 

(«i^J 



(«i^J 



NL) 



4 



region 1 
region 2 
region 3 
region 4 

region 5 
region 6 

region 7 

region 8 



(44) 



In we have used the standard "smectic penetration 
length" A = \/ Kq/Bq, where Kq and Bq are the "bare" 
smectic layer bending modulus and compression modulus 
respectively; i.e., the bending modulus and compression 
modulus for the smectic in the absence of disorder. The 
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loci separating the regions in Fig. [4] are: 



FH: q.^^L^iq^i^L? " (45) 




where go = 2^0, = (C^^j'/A, A = C/(l + W2), 
and a is the smectic layer spacing. Here the positive 
definite universal exponent 773 is defined in Eq. (jB3| . We 
have calculated it in the e = 5 — d-expansion and found: 

^3 = l + 0(,2). (52) 
5 
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FIG. 4: Illustration of the eight regions listed in (144 |l . each 
of which exhibits different scaling of the light scattering in- 
tensity in the A phase as T — * T\q. The light scattering 
intensity is dominated by random field fluctuations in regions 
f and 2, by random tilt fluctuations in regions 3, 4, 6 and 7, 
and by random compression fluctuations in regions 5 and 8. 
Above locus OFH the common denominator in formula (|42p 
is dominated by Bq^, while below OFH it is dominated by 
Kq\ in region 6 and by -Dg^ in regions 1 and 3, respectively. 



Of course, none of these crossover lines is sharp, since 
none of the crossovers is abrupt. Rather, as is the nature 
of all crossovers, these are loci near which two terms start 
to become comparable. Hence, these crossover bound- 
aries are only defined up to 0(1) factors. 

The fairly complicated scaling behavior embodied in 
(H^ simplifies considerably when we take gz = and 
vary q± (i.e., move along locus 1 in Fig. [4]). The sec- 
ond term in (|42| then vanishes. For g_L ^ , Kq]^ 
dominates in the denominator of the first term in (|42p. 
which term itself dominates the third term, and using 
dMl) we obtain {N^ {q) Nj- {-q)) oc g-4+2'?K -r,* ^ 

qi_ <C ^2^, the first term dominates at short wavelength, 
and the third term dominates at long wavelengths. The 
crossover between these two is at 



= (efr ' oc (53) 



is much larger than as T ^^ac' provided 
2> + r]t — > 1, a condition we can show to be sat- 
isfied. Clearly, for q ^ (^j^)^^, if there were no anoma- 
lous elasticity of K (i.e., if rjK — 0), the light scattering 
lineshape would be a Lorentzian squared, with correla- 
tion length In reality, the lineshape is quantitatively 
different (e.g., its tails scale like qj*^"^^"^^^*), but still 
qualitatively like a Lorentzian squared, with linewidth 
^2^. For q <C (^;5^)~^the light scattering diverges as a 
power law. More precisely. 



[ [D{T) + K {q^,q, = 0, T) qlf [D{T) + K {q^,q, = 0, T) ql] q^ql ) 
< H^^l) {^) , (ef) «q^«C (54) 

_L / ± \ 
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Comparison of this expression with hght scattering data 
should allow easy determination of (T) and the combi- 
nation of exponents 2r]K — Vt- Fitting the T-dependence 
of U{T) to (T - TAcy^ then determines z/j.. 



Another simple locus is obtained by taking q± fixed 
in the range [Tl[ qj_ ^ S'Hd varying q^ (i.e., 

moving along locus 2 in Fig. |4]), which gives 



9. » ('Zie^L)^ (e^i)"' 



(55) 



7^ ^ j g± < < (6 



< < 



where qf^ ~ (Ctvl/'C-l) ^ • From this locus 
can be read from the point where a log-log plot of / (q) vs 
qz for fixed changes slope from — —4 to —4. Again, 
once £,z(T) is known, can be determined. Note that 
measurements of the exponents ^ -I- 4 and ^ -f 2, taken 
in conjunction with knowledge of rjt — 2rjx from locus 
1, determine rjK, T]t and rjc (recall that the anisotropy 
exponent ( = 2 - ^). 

To summarize our light scattering predictions for the 
v4-side of the critical regime: measurements along locus 
1 and 2 can be used to determine the exponents i'±, i>z, 
rjK, and 774, as well as testing that rjB = (which was 
assumed in the above discussion). These exponents can 
then be quantitatively compared with our predictions, 



and, furthermore, test the exact scaling relation (|41|). 
and combined with specific heat measurements to test 
the exact scaling relation ([39]) for a, and with the afore- 
described dielectric measurement determination of (3 to 
test the equally exact scaling relation (j40|) . 

We now turn to the smectic C (T < T^c)- For q±^± » 
1 or qzS^z ^ 1, the light scattering is the same as in the 
high temperature case, since this range of wavevectors 
corresponds to the critical regime, in which one cannot 
tell whether one is on the high- or low-temperature side 
of the transition. In the other limit - namely, when both 
q±^± <C 1 and qz^z ^ 1, we are in the smectic C-phase 
range of wavevectors, and Cij is given by: 



/ 



A,, {q',T)qlql 



A.' {T)qlql 



\ 



\ (7(9 T)ql, + Bql + K {q', T) q^ 

A^,iq\T)ql,ql 
{j{q'',T)ql,+Bql+K{q',T)qt, 



7(9', T)ql, + Bql + K iq',T) qj, j J 



(56) 



where qx' , qs', and qz' are defined in Eqs. ([HI El flUl) 
respectively, 7 and Agi^^' are anomalous and given 
in Eqs. ([UlTlin]), and Az' is wavevector g- independent. 
In addition, the temperature dependent quantity r(r) 
defined in the coordinate transformations Eqs. ([51[^[TU|) 
has the critical scaling: 



r(r)(x 



NL 



(57) 



Since the azimuthal symmetry about the z-axis has been 
broken due to the tilting of the layers, the light scat- 



tering also loses azimuthal symmetry, and becomes fully 
three-dimensional. It is most conveniently described in 
the transformed, non-orthogonal wavevector coordinates 
(i.e., (f'-space), whose relation to the lab wavevector co- 
ordinates (i.e., (f-space) are given in Eqs. ([51 [51 [TU)) . 

As in the A-phase, in the C-phase the light scattering 
displays many different regimes of very different scaling 
with wavevector and temperature, due, as there, to both 
the crossovers between domination by different types of 
disorder, and by the crossovers between dependences on 
different components of q. Here, the three-dimensional 
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q' space is divided into the five regions illustrated in Fig. [5l in which the light scattering intensity is given by 
I 



A 



^KihL?{-it:T '""^'(^^9.')^""'^^'^^'^/'^' (^)', region O-EDIG 



region 0-ABCH 
region 0-AFEC 
region 0-CEGB 
region 0-FJDE 



(58) 



where we've defined the C-phase crossover exponents: 











H 




B 


G 
















6^; 






1 



FIG. 5: 3D picture in tlie transformed reciprocal space (g" - 
space) illustrating the five distinct regions with different 
wavevector-dependences of Cij{q) in the C phase, which is 
given in ([58}. 



The crossovers between these five regions obey 
OAC : q,,^, = 



A 



OFE : q,ji, = {qs'U)^'' 



'73/2 



(59) 
(60) 




OGE : q,,(, - 
ODE : - iqs'C^)^'-' 

OEC : fe. = (g.'6)'^='^' 



(61) 
(62) 
(63) 
(64) 



(65) 
(66) 
(67) 



and the quantity A° = Aj! + A^^^ {^nl) /wqKo, where 
go and are the "complete bare" values of of the coeffi- 
cients of the anharmonic terms in the Hamiltonian (1831) . 



Now let us narrow our discussion to wavevectors in the 
Qz-Qx plane and measure Cij{q) given by 



- (^-;-^.^' + A.,(,-r),go,g ^ ^^^^ 



jiq,T)ql, + Bql 



where we alert the reader to the fact that this expression 
contains both the rotated wavevector component qx' and 
the un-rotated wavevector component qx- This is not a 
typo! 

This two-dimensional qz-qx plane can also be di- 
vided into different regions with different wavevector- 
dependences of Cxx{q), including (for illustration) the 
critical regime (i.e., qz ^ or q± ^J^)- (Note that 
we are now working with the un-transformed q' s). This 
leads to the 10 regions which are illustrated in Fig. [6l In 
these regions, the light scattering results are: 
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Cxxiq) 



A 



region 1 

region 2, 3 

region 4 

regions, 6 

region 7, 8 
region 9, 10 



(69) 



At short wavelengths (i.e., Qx ^ Cl^ or ^ Cj^); that 
is, in the critical regime, the crossovers between the dif- 
ferent regions remain the same as those of high temper- 
ature phase, which were given in equations (|45l [46|) . At 
long wavelengths (i.e., both ^ and qx ^ ^J^) the 
crossovers are given by: 



FIG. 6: Illustration of the 10 distinct regions in qz-Qx plane 
with different wavevector-dependences of Cxx{q), which is 
given in (I69|l for the C phase. Line OA is define by q^i = 0. 



OE, OE' : = 

OF, OF' : (9,6)^^' = (qxU)'^" 



||(|)(g.-rg.) 



Two possible experimental loci through the qx-qz plane 
are shown in Fig. [5) By holding 3, qx fixed and varying 
qz (locus 3), we can determine — fjx' /Cz' — 4 from a log- log 
plot of I{q) vs qz, and also deduce ^z from the location 
of the point at which the log-log plot changes slope from 
-fjx'/Cz' - 4 to -T]t/C - 4 . Likewise - {rix' + 2fi^) /C,x' 
and can be determined by moving along locus 4 (the 



Hard Direction 




To completely determine the anomalous exponents for 
the low temperature phase, one more experiment is nec- 
essary. One possibility is to measure the light scattering 



at both (/z = and qx — Q, varying qa, for which (j56[) 
leads to 
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Css iq) 



A 



?7t-2r;jf 



r 



(70) 



Such measurements along this experimental locus deter- 
mine 277x — 77s' — 4 and also measure Combining this 
result with the data for —fj^i/Q/ — 4, — {fix' + '2f]-y) /Cx' 
and the scaling relations (US]), we can deduce the expo- 
nents ffj, f/K, fjs' 1 fix' and, furthermore, check the scaling 
relation 1151). 



The remainder of this paper is organized as follows. 
In section |lTl we derive our model for the smectic in 
stretched aerogel, and perform the "replica trick" on it. 
In section Unl we derive the renormalization group recur- 
sion relations for the model near the AC transition, iden- 
tify their stable fixed point, and calculate the thermody- 
namic critical exponents in a d = 5 — e expansion. In sec- 
tion |TVl we calculate the critical exponents. In sections 
IVl and IVII|. we calculate the wavevector-dependences of 
the elastic constants and disorder variances. In section 
IVIl we treat the smectic A and C phases of the model 
and show that the former is in the "XF Bragg glass" 
universality class, while the latter is in the "m = 1 smec- 
tic Bragg glass" 's'l universality class. In section IVIIII 
we calculate correlation functions near critical region. In 
section HXl we discuss the validity of our predictions for 
the critical behavior against the unbinding of topological 
defects and the loss of orientational order. 



II. MODEL 

As in earlier treatments of smectics A and C[T6| in 
clean environments, we expect, even in the presence 
of disorder, that the important fluctuating variables in 
our system the "layer displacement field" u(r), and the 
molecular axis fi{r). The layer displacement field u{r) 
is defined, as usual, as the displacement along the mean 
layer normal of the layers relative to some ideal refer- 
ence configuration of flat, uniformly spaced layers. We 
will choose the orientation of these reference layers to be 
perpendicular to the stretch axis. 

We will derive our starting Hamiltonian by beginning 
with that for a smectic A in a clean environment, and 
then modifying it to reflect the effects of the aerogel, and 
to allow the molecular axes h{r) to tilt away from the 
layer normal. 

The Hamiltonian for the smectic A phase in a clean 



environment is |16l |: 
'K 



H, 





4" 







f (Vi«)^ + |(5,zi)2-|(5.^)|V^ur 



(71) 



where, in the second equality, g — w — B. This equality 
of w, and B, which is, of course, a consequence of sim- 
ply expanding out the square of the B-term in Eq. (|71|). 
is also a consequence of the global rotation invariance of 
a smectic in the clean system, since it is that invariance 
which forces the B-term to be precisely the square of the 
rotation invariant quantity E = dzU — ^|V_LUp in the 
first place. 

Once the smectic is placed in a non-rotation invari- 
ant environment like stretched aerogel, we no longer have 
the constraint of rotation invariance. However, obviously, 
any term that was allowed by symmetry in the rotation 
invariant case will certainly be allowed when the sym- 
metry is lowered. Hence, the B, g, and w terms in the 
second equality of Eq. ([71]) are still allowed when the 
smectic is put in anisotropic aerogel. But since that sys- 
tem is no longer rotation invariant, there is no longer any 
symmetry requiring that B = g = w. Thus, for a smec- 
tic in anisotropic aerogel, part of the Hamiltonian will be 
just the expression after the second equality in ([71]) with 
no constraints relating B, g, and w; that is, all three will 
now be independent parameters of our model. 

There are additional terms not present in (|71|) that 
are allowed in our Hamiltonian due to the absence of 
rotation invariance in our problem. In particular, there 
can, and hence, in general, will be terms picking out a 
preferred orientation for both the layer normal N and 
the molecular director n. We will assume that both of 
these preferred directions are along the stretch axis z. 

The simplest coupling that will accomplish this is: 



H, 



stretch 



;^2 



d^r M{N ■ zy +Q{n- z) 



(72) 



with M and N both > 0. 

In addition, the molecular axis fi and the layer normal 
N should couple. As noted by deGennes[lB|, the simplest 
such coupling takes the form: 



H 



Nn 



P / d''r\N - fiY 



(73) 
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Using the well-known [ly| geometrical relation 



N : 



(74) 



between TV and the layer displacement field u, defin- 
ing 5n = n~ £, and dropping terms of 0(|V±m|'*, |(5n_L|^) 
(which can be shown to be absorbable into a shift of the 
quartic coupling w in Eq. (|7ip ). and adding terms reflect- 
ing the randomness of the aerogel, we obtain our starting 
Hamiltonian: H = Hi + H2, with: 



+ - v.. 

and 

Ho = / d^r 



+ h{r)-Vu + Vp{u-(t){r)) 



(75) 



M\V±_u\^ + P\V^u + + Q\5rn_\ 

(76) 

The h[r) in (j75p is a quenched random field which we 
take to be Gaussian, zero mean, and characterized by 
short-ranged anisotropic correlations 



h^ [r] h, (f ') [AtS^^ + A,<5f,] [r - r') . 



(77) 



One might reasonably question this assumption of short- 
ranged correlations of the disorder, given the known frac- 
tal structure of aerogel. Indeed, while reference |13j] ar- 
gued that for liquid crystals in aerogel the components 
of h orthogonal to z do not have such long-ranged cor- 
relations, hz{r) does. Hence, our results are not directly 
applicable to these systems, although they could be ap- 
plied to systems of higher aerogel density, which would 
have shorter ranged correlations. 

The field (/<(r) in ([75)1 is also a quenched random field 
with only short-ranged correlations, and is uniformly dis- 
tributed between and a, the smectic layer spacing. The 
function Vp{u — (f>) is periodic with period a. 

The physical interpretation of the quenched random 
fields h{r) and Vp{u — (j)) is very simple. The random 
field h incorporates random torques and random com- 
pressions, coming from the _L and z components of h, 
respectively. The function Vp{u — <p{r)) represents the 
tendency of the aerogel to pin the smectic layers in a set 
of random positions (/'(r), modulo the smectic layer spac- 
ing a, which is why Vp is periodic in its argument. The 
Hamiltonian ([75| is identical to the elastic theory of a 
smectic A in isotropic aerogel developed in , with one 
crucial exception: in a smectic in isotropic aerogel, rota- 
tion invariance requires that g = w = B, while for our 
problem, even at T = Taci where softness is recovered, 
g and w are still free [9], because rotational invariance 
is still broken. The remaining effects of the anisotropy 
are incorporated in the M and Q terms in Eq. ([7S|) . It 
is the P term that actually drives the AC transition in 



this model. When P is positive, it is energetically fa- 
vorable for n to be normal to the layers (which implies 
6n = —W±u), while, as we shall see in a moment, when 
P is sufficiently negative, it is energetically favorable for 
the molecules to tip relative to the layers, thereby putting 
the system in the C phase. 

To see this, note that after a linear change of variable: 



Sn' ± = Snj_ + R\7 ±u 
with R = P/{Q + P), (I7HD becomes 



H2= d" 



D{T) 



^±u\' + {P + Q)\5n'^\^ 



(78) 



(79) 



with D{T) = 2M + 2QP/{Q + P). Assuming when tem- 
perature drops both M and Q remain positive and P 
decreases continuously from positive to negative, which 
drives the system from A phase into C phase, D — 
M + QP/{Q + P) becomes negative first. Therefore 
just below Tag, at the new ground state, |Vj^it| is 
non-zero, Sn' ±^ — 0, which, combined with (|78p . gives 
dn± = — _RV_lu. This implies that both the layers and 
the director tilt, but in opposite directions. Clearly, if 
Q ^ M, the layers tilt more than the molecules, while 
in the opposite limit {M ^ Q), the director tilts more, 
as observed in [l^ . 

Since the fluctuations of Sn' ± are massive, we can set 
dn'± = in and use the sum of (jTS]) and (HH), which 
now only involves u, as our complete Hamiltonian: 



H 



d'^r 
((92u)|V^u 



+ h{r}-Vu + Vpiu-(l){r)) 



(80) 



Since |V^m| is proportional to the tilt angle of the 
smectic layers, the coefficient D{T) is positive in the 
A phase (favoring alignment of the layer normal with 
the stretch axis), and negative in the C phase (favoring 
tilt of the layers). Hence, by continuity, at T = Tac 
D{T) vanishes. In what follows, we will assume that 
D{T) (xT - Tac near Tac- 

We can now treat the quenched disorder in this Hamil- 
tonian using the replica trick [l3| . This trick starts 
by assuming that the actual free energy and correlation 
functions measured in an experiment on a single sam- 
ple will be close to the average values of these quantities 
when averaged over many realizations of the disorder, 
with a weight for this average given by the distributions 
described earlier. To compute, e.g., the average Free en- 
ergy F — — Tlog Z, where Z is the partition function, we 
use on the identity logZ = lim . can now be 

n^O " 

computed by doing a repeated functional integral over n 
"replicas" of the field u: 
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Dui e-'3^^("i.'^.'A) / Du2 e-P^i-^^'^- 



Dun e-Z^^t""'''- 



J J Dua e 



(81) 



a=l 



The advantage gained via the rephca trick is that the 



average over the random fields h and (j) (which, we note, 
are the same for all n replicas), can now be done first, 
leading to an effective "replicated" Hamiltonian for the 
set of n Uc's alone. The average over the random fields 
h and (j) couples the previously uncoupled replicas. 



Details of how this averaging process is performed can 
be obtained in reference [l3[. The calculations here are 
virtually identical, the only differences being the slightly 
different form of the starting Hamiltonian (1751 176]) . After 
replicating and integrating over the disorder h{f) utiliz- 
ing Eq. ([77|l we obtain 



J 



H[Ua] 



2 , W 



K (Viii„) + B {d.Uc^Y - g{dzUo) (V±wa) + - 



D{T)\V^u^ 



E 

a,/3=l 



(82) 



where Ap(mq, — up) is a periodic function with period a, 
the smectic layer spacing. 

III. RENORMALIZATION GROUP 

In this section we derive the RG flow equations for 
studying the phase transition. The basic idea of the 
RG is to eliminate the short-length degrees of freedom 
of u and see how this affects the long-length physics. At 
length scales shorter than the correlation length, since 
the tilt term in Hamiltonian ([82)1 is negligible compared 
to the bend term, our model is effectively similar to that 
for smectic A in isotropic aerogel, in which the tilt term 
is absent due to the symmetry. In the latter, reference 
showed that the random field disorder /S.p{ua — up) is 



irrelevant compared to the random tilt disorder in d < 5. 
Repeating the (virtually identical) calculation for this 
problem, we also find that /S.p(ua — up) is irrelevant in 
our problem at length scales shorter than the correlation 
length. However, it becomes relevant at longer length 
scales in the A phase, which will be discussed in sec- 
tion IVII In addition, the random compression term Ac 
can also be shown to be irrelevant in the renormalization 
group sense. However, while it does not affect the renor- 
malization group recursion relations to lowest order in e, 
it does prove to be important for the correlation func- 
tions. We will therefore treat it in detail later in section 

Eul 

With these simplifications , we can study the AC phase 
transition using the following truncated Hamiltonian: 



A r " 



K (Viua)' + B {d.uc.)^ - g{d,Uo.) (V^Ua)" + j 



'V±Ua 



DiT)\V±Ua\ 



Up. 



(83) 



a,/3=l 



This Hamiltonian's nonintcracting propagator Gap{q) =< Ua{q)up{—q) >o is found to be 



G^p{q) = TG{q)S^p + AtqlG{q)^ 



(84) 
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modulus B 




FIG. 7: The graphical correction to B coming from the com- 
bination of the two cubic vertices. 



with 



G(q) 



(85) 



To qualitatively estimate the effect of the anharmonic 
terms in Hamiltonian (j83p . we perform the ordinary per- 
turbation theory. We find the perturbation theory breaks 
down in d < 5. This can be seen by calculating the one- 
loop graphical correction to B, which is represented by 
the Feynman diagram in Fig. [71 An analysis of the Feyn- 
man diagram gives 



_ 9_ 

AT 

in which 



2 n n „ 

EEE'?'""('?>M-9) (86) 

a=l/3=l q "'P 



-f A2piG(p)4 , (87) 

where 5af3 is a Kronecker delta function. This diagram 
leads to two contributions. One of them is 



9^ 
4T 



n n „ 

T.T.Y.'l'''-('l)M-Q) p1G{p)\ (88) 

a=l/3=l ''P 



which is actually a contribution to the random compres- 
sion. We will not discuss this contribution here since the 
random compression is irrelevant. Later we will come 
back to it when we calculate the wavevector-dependence 
of Ac- The other contribution is 



a— 1 q 



X / pi [TG{p)^ + 2AtplG{m , (89) 
Jp 

from which we obtain the correction to the compression 



6B 



p1 



],{Bpl+Kpir 
dp. d''-'p± 



2tt Ji. (27r)'i-i {Bpl + Kp^ 



3 Sd-iBAt /g\^fB 



{2nf-' {5 ~ d) ^BJ \K5 



1/2 



L 



5-d 



(90) 



where, in the second equality we throw out the thermal 
fluctuation piece because it is less divergent, Sd-i is the 
surface area of a d — 1 dimensional sphere with unit ra- 
dius. This result diverges as a power law of the system 
size L right at the critical point in d < 5, which implies 
that the upper critical dimension due for this phase tran- 
sition is 5, below which Gaussian theory breaks down. 
Since the physical dimension d = 3 is less than due, the 
standard momentum shell renormalization group (RG) 
transformation has to be employed to study the phase 
transition. 

Because a few features of our treatment are non- 
standard, we will describe our approach in some detail. 

The flrst unusual feature is that we use an infinite 
hyper-cylindrical Brillouin zone: \q\ <A, —00 < < 00, 
where A ~ 1/a is an ultra-violet (UV) cutoff. We sepa- 
rate the displacement field into high and low wavevector 
components. Ua (r) = u'^ (r) + {r) where m„ (r) has 
support in the hyper-cylindrical shell Ae~^ < q± < 
—00 < qz < 00, and (f) has support in the remainder 
of the hyper-cylinder (i.e., q± < Ae^^, —00 < q. < 00). 
We then integrate the high wavevector part m> (r), and 
rescale the length and long wavelength part of the fields 
with f± = r^e^, z = z'e'^^, and (r) = e^^u'^ (r ') so as 
to restore the UV cutoff back to A. This rescaling leads 
to the zeroth order RG ffows of the effective couplings 



K{i) 

B{£) 

At 
T 



At 



T 



(d+Lj+2x~3)e 



From these dimensional couplings one can construct four 
dimensionless couplings: 



91 = 




-1 {B/K' 


)-'A'^-^, (91) 


92 = 


Cd- 


-1 (B/K^ 


)^AtA'^-5, (92) 


53 = 


Cd- 


-i{9/Bf- 


(B/K^)^ AtA'^-^ (93) 


94 = 


Gd- 


-,{w/B){ 


B/A'S)^ AtA''-^ (94) 


i-l = 


Sd- 




\ The RG fiows of the di- 



mensionless couplings are independent of the arbitrary 



FIG. 8: Graphical corrections to the bend, tilt and random 

tilt terms. FIG. 10: The graphical correction to the linear term 2, 9z 




FIG. 9: A graphical correction to the tilt term. 

rescaling exponents and given, ignoring graphical correc- 
tions, by 

51 W = gie^^-"'", (95) 

52 W = .92e(5-'^)^ (gg) 

53 W = 53e(^-'*)^ (97) 
54(^) = gie^'-'^'. (98) 

gi is just the diniensionless coupling found by Grinstein 
and Pelcovits[9] in the pure smectic problem. It becomes 
marginal and leads to very weak anomalous elasticity in 
d — 3. 52 is actually the dimensionless coupling found by 
Toner and Radsihovsky Ts^ in the isotropic disordered 
smectic problem, where it becomes relevant in d < 5 and 
leads to much stronger anomalous elasticity in d = 3. 53 
and 54, the truly nontrivial dimensionless couplings in 
our problem, also become relevant in d < 5. As we will 
see below, they control the structure of the RG flows and 
hence the critical behavior. 



FIG. 11: A graphical correction to the cubic term. 




FIG. 12: Another graphical correction to the cubic term. 
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FIG. 15: Another graphical correction to the quaxtic term. 



FIG. 18: One more graphical correction to the quartic term. 
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FIG. 19: The final graphical correction to the quartic term. 



(i.e., ii'-term), the tilt term (i.e., D-term) and the ran- 
dom tilt term (i.e, A^-term). The Feynman diagram for 
these corrections are presented in Fig. [8]and[9l Following 
Feynman rules, from the diagram in Fig. [5] we obtain 

2 n n 

a — 1 f3—l q ij 

X / pi {Pf + it )pI Ga(3 {p)Gapip + q}, 
Jp 

(101) 

2 n n 

ct — 1 /3— 1 q ij 

X r pipfplGo.p{p)Go.p{p + q), (102) 
Jp 



Since the integration over the high wavevector com- 
ponents of Ua is always performed in the momentum 
shell (— oo < Qz < oo, Ae^'^^ < q± < A), it can be ac- 
complished perturbatively in nonlinearities of H[u\. The 
change of the effective Hamiltonian of the long length 
modes due to the integration is 



SH[i 



{Ha[ 



^ / tt3 



6T 
1 



24r3 



(Hi 



(99) 



where GafS is given in ([57)1 . Note that we have kept the 
(f-dependence of the internal legs (i.e., the g-dependence 
of the integrands in (jlOip . (|102p ) since it is necessary 
for obtaining the contributions to the bend term. Taylor 
expanding the two integrands around q — 0, we find that, 
0(1) terms lead to the corrections to At and D; 0{q±) 
terms vanish after the integration; 0{q^) terms lead to 
the corrections to K. 

Let us first calculate the corrections to A( and D. Set- 
ting q inside the two integrands to 0, (jlOip and (|102p 
become identical and can be added together to get 



where Ha represents the anharmonic terms which in- 
cludes the cubic and quartic terms, the average over the 
short length modes is done by only using the harmonic 
part of the Hamiltonian. 

Now we will discuss 5H[u'^] to lowest order (one loop 
graphical correction) in detail. 

Since as shown in ([M]) . the calculation of SH[u^] is just 
a modified perturbation theory, the evaluation of the one- 
loop graphical correction to B for long wavelength modes 
can be accomplished by modifying (j90p . restricting the 
integral range within the momentum shell {—oo < q^ < 



oo, Ae 
SB = 



-de 



< q± < A). Doing so we obtain 
2 r> 



^ I pl[TGip)-' + 2A,plGim 



Pa 



{Bpl 
dpz 



2-K 



Kp\Y 



PA. 



pI 



16 



Cd-iB/^t 



2 



(27r)'i-i [Bpl -f Kp\Y 



b) 



1/2 



--Bg^dL 
16 



(100) 



S'^ 5Z 51 XI XI 1ilf^o,{q)u(s{-q) 

a=l/3=l q ij 

X r pip'rplGo.p{p)G^p{p), (103) 



from which we obtain two distinct contributions. One of 
them is the contribution to the random tilt term with the 
corresponding correction to the random tilt variable At: 



5At « '-^ j\lplAjG{p)' 



,2 r> 
.2 



-Ai 



-A 



pIpI 



{Bpl + Kp\Y 



dpz 
2n 



P± 



p\pI 



(27r)rf-i {Bpl+Kpl 



4 ^3 



(104) 



Now we discuss graphical corrections to the bend term 



The other one is the correction to the ordinary tilt term 
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with the corresponding correction to the tilt modulus D: 

»> 



{5D), 



V 



A 2 



^ 16 32 '^^ ' 



(105) 



where we have used {SD)i instead of {6D) since there are 
more corrections to D from other Feynman diagraphs, 
which will be calculated later in this section. 

Now let us calculate the corrections to the bend term 
from (|10ip and (|102p . As mentioned earlier, we have to 
Taylor expand both of the integrands and use the 0{qj^) 
term. Fortunately, since p02p is part of poip . we only 
need to do the calculation once for (|10ip . which, however, 
is still very involved. In the following we will describe the 
calculation in detail. 



J 



Only keeping the most divergent pieces which carry 
the factor Sap , the integrand in (|10ip is simplified as 



pi (pi + (ii)pl 



AtP\5ap 



{Kp\ + BpmK\p + cr\^ + Bpl) 



{Kp\ + Bpl){K\p + <f\^ + Bplf 



(106) 



(107) 



For convenience, we rewrite K\p-\- c^'^ as 
K\p + q\^ = K{pl + f{q^)) 
with ,f{qi_) defined as 

/(gl) = 2pl{2p^ ■ 9l + ql) + {2p^ ■ qi_ + qlf . (108) 
Now let us do the following two Taylor expansions: 



1 



1 



{K\p^ + q^\^ + Bpl 



Bpl 



1 



Kp\ 
1 



(K|p_L + + Bp; 



2\2 



{Bpl+Kp\f 



Kf{pi_) K^.f{p±) 
Bpl + Kp\ [Bpl + KplY 

2Kf{p^) 'iK^f{p 
Bpl + Kp\ {Bpl 



1 - 



Kp\f 



r 



Plugging these two expressions into ()106p makes it ex- 
tremely long and complicated. For better management. 



We divided it into four pieces: 



pipjpl^tp\5aii 

^2^3 



{Kp\ + Bplf 

>lqj-plAtp\Sa/3 



,2« 



(Kp^ + Bp 
pipfpl^t\p + q\^Sap 

(Kpi+Bpir 

Phj-pl^t\p + q\^Sal3 



{Kp\ + Bplf 



Kf{qi_) 




Bpl + Kp\ 


{Bpl+Kp^^)^_ 


Kf{q±) 


K'Piq^) 1 


Bpl + Kp\ 


[Bpl+Kp^^rl 


2Kfiqj_) 


, iK'nq±) 1 


Bpl + Kp\ 


^ {Bpl+Kpi)^\ 


2Kf{q^) 


, 3K'fiq±) 1 


Bpl + Kp\ 


{Bpl + Kp^^y\ 



(109) 
(110) 
(111) 
(112) 



Inserting Eq. (|108p into the above four pieces, we obtain O(g^) terms: 
I 



hiqi 



-2KAtql5ap 



2 4 1 I 

PzPlprPj 



[BP^+KpX)^ 



^plp\pipfPkPi 



k I 



[BP] 



KpiY 



-lQK^At5apY.Y. 



^plplpipjpjpi 



{BP^ + Kp\f 



(113) 
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from p09l) : 



from PTU)1 : 



from (fTTTjl : 



= -4KA,q^6^,J2l^j0^^t^^ (114) 



2 4 11 ? ? I I I 

PzPlPtPj ar.-^ A ^^^^PzPiPtPjPkPe 

,4 \4 



from (|112p . Finally the total correction to the bend term 
can be written as 



2 n " r> 

^ E E E E Qht^c.(<i)M-'i) X y /fe(gi) , 

(117) 



a— 1 /3— 1 g ijk 



where the subscript k — 1,2, 3, 4. 

The main trick involving doing the integrals in (|117p 
can be summarized by the following two integrals: 



(118) 



pW.np') = :jj6,, / ip'mp') 



pWkPW) 



1 



p' 



d'{d' + 1) 



(SijSke + SikSji + SiiSjk) 



X / (y)^r(p'), 



(119) 



where both function T{p') and the integral region are 
spherical symmetric, d' is the dimension of the space in 
which p' lies in. 

Now we plug in the expression of I/, into (|117p and do 
the integrals. Let us start with Ii, which has three terms. 
Integrating the first term gives 



9 4 11 
PzPlPiPj 



(i?P2 + Kpi)^ 
2KA, f< pIpI 

^_^q.Oo,p0^j I ^Rp2_^^„4^,4 



p [BP^ + KpiY 



pIp\ 



— -q±do,p5,j / dp^ / dp±- , 4 ^4 

d-1 Jhe-it {BPi+Kp\Y 



(120) 



integrating the second term gives 



-4i^Ai<5„^EE9fc9^ / 

1. I) J V 



2^2 



i_ i_ / PzPlPiPjPkPe 



k i 



AKA 



^2 _i ^o'P E E^'^*-'''^'''^ + SikSji + 5u5jk)qk(li 
k e 

(BP^ + Kp^^)* 
X E E^'^*^''^''^ + '^ife'^i^ + Sii6jk)qkqiCd-iA^~'''d£ 

k £ 

3 

-Y^A, (^-^) ' <5„^(<5.,gi + 2q^q^)Cd-iA'-''di; 



integrating the third term gives 
16K'AtS^f,J2J2lkll I 



k i 



p\p\pip\PkPi 

[BPl + Kp\f 



k I 



Pip\0 



p [BPj + Kp 



4 u 



5 / 1 \ ' 
X E '^{^■tj^M + kk5jt + 6ti5jk)qtqiCd-ih- 

k I 



334 At ) 5^p{5,,qi + 2g,^g/)Cd-iA^-'^c;^ . 
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Adding the three results together we get 



hiqi 



\BK 



{121) 



Plugging this result into (|117p we find corrections to the 
bend term: 



the corrections to D are illustrated in Fig. [8] and [9l re- 
spectively. The former has been discussed previously and 
leads to {5D)i given in (jlOSp . A simple analysis of the 
latter gives 



p1 



3 r> 

4"" {Bpi + Kpl + Dpl)^ 

(^^KA^ ^ —D)g^d£. 
46 32 



(126) 



2 n n 



1=1/3=1 q ij 



1 

64^ 



256 \BK J 

a q 



(122) 



from which the corresponding correction to the bend 
modulus K is 



'l28^* ^Bk) 
128 ^ 



(123) 



Naively one would treat the sum of {5D)i and {SD)2 as 
the total correction to D; however, a more careful analy- 
sis shows that this treatment is incorrect. One has to take 
into account that the disorder averaged ground state of 
the rescaled system has been renormalized under the RG. 
The Feynman diagram in Fig. 1101 shows a contribution 



(127) 



a=l 



which is linear and does not exist in the original Hamil- 
tonian. This implies that the disorder averaged min- 
imum of the Hamiltonian has shifted from dzUa = 
to dzUa — —6L/B. Therefore, after each cycle of the 
RG, the Hamiltonian needs to be reexpanded around 
the new minimum, which is accomplished by changing 
variable, defining u'^ = Ua + {5L/B)z. This treatment 
then makes the cubic term —g{dzu){y ±uY /2 to create 
another contribution to the tilt term {5D)^\V ^u']"^ /2, 
where {5D)^ = gSL/B. 6L can be easily calculated from 
the Feynman diagram in Fig. 1101 



Similarly, inserting l2{<li), hi^x) and hiq^) into pi7p . 
we obtain three more corrections to K: 



1 

32 



Kgsde, 



m)3 = Y^I<93di, 



(124) 



Now we have found all the correction to the bend modu- 
lus K from (HnH). 

As for (|102p . since it is part of (llOip . its correction to 
K can be easily shown to be {SK)i + {SK)^. Thus the 
total correction to the bend modulus K is given by 

SK = {5K)i X 2 + {6K)2 + {6K)3 x 2 + {SK)^ 



= -^^Kg,di. 



(125) 



Now let us turn to the calculation of the graphical cor- 
rections to D. The two Feynman diagrams which give 



6L 



1 



p1 



p 



[Bpl+Kp\ + Dp\Y 



« [-lKk^ + ^D)[^^g,di. (128) 

Taking {SD)^ into account, the total correction to D 
should really be 

5D = {5D)i + {5D)2 + {5D)z 

= ^(53-<?4) + ^if(,g4 - 53). (129) 

In the isotropic disordered smectic A problem Q the rota- 
tion invariance prohibits the tilt term and this symmetry 
should be preserved under the RG, which provides an ex- 
tra check of our calculation. Imposing 53 = g^ (restoring 
the rotation invariance) in (|129p . we find indeed that 5D 
vanishes. 

Now we calculate graphical corrections to anharmonic 
terms in the Hamiltonian ([55]) . This calculation has been 
avoided in the isotropic disordered smectic A problemQ 
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by using the symmetry argument that the global rotation 
invariance requires 

SB = Sg = Sw . (130) 

This argument does not apply in our problem due to the 
symmetry breaking; however, it can be used to check our 
calculation. By imposing B ^ g = w, our results should 
also satisfy (jl30|) . 

All the one-loop Feynman diagrams which give the cor- 



rections to the cubic and quartic terms are shown through 
Fig. [T4]to[T9l The corresponding graphical corrections to 
g and lu are in order referred as {Sg)i, {Sg)2, iSg)^, ((5w)i, 
{Slu)2, {610)3, {6^)4, (<^^)5, {Suj)6i respectively. Since the 
evaluations of these diagrams are all similar, we choose 
to only show in detail the evaluation of the Feynman di- 
agram in Fig. [191 which is the most complicated, and for 
the rest we only give results. An analysis of the Feynman 
diagram in Fig. [19] gives 



T3 V 



00 00 00 00 



EE EE E E '^(«^i + '?2 + g;3 + g3)(i9n)(ig2'j-)(*'?3'i)(*'?il) 



a=l 13=1 7=1 ri=l <fi ,13*2 ,93 ,g4 



(131) 



r 



in which 



^ AtqlG{q)'){TG{q)60^ + AtqlG{q)') 

AtqlG{q)^){TG{q)So.rj + AtqlG{q)'). 

(132) 



{TG{q)S^p . 
{TGiq)Sy^ - 



Note that the quartic term in the Hamiltonian (|83|) only 
involves one sum over the replica index, while ()13ip has 
four. Thus only those terms in GapGpjGjriGari which 
carry at least three Kronecker delta functions can lead 
to the contributions to the quartic term. Plugging (|132p 
into (|13ip and keeping the most relevant pieces, we obtain 



iSLo)e = -Sg'At / piptG{q)' 
Jp 

32 gi 



(133) 



The corrections to g and w from other Feynman diagrams 
can be calculated in a similar way and are 



(-5.9)1 
(-55)3 



{Sg)2 = -—g^gdl, 



16 



g4Ujd£ , {Slu)3 



3 

16 
3 

32 



g^ujdi , 



16 



93^d£ , {Suj)5 



16 



gsusdi . 



Adding them together we get 
3 9 



Slo 



3 gl 

—L 

32 54 



15 
32 



giUj 



rff3W. 



(134) 
(135) 



Performing the rescalings described earlier and using 
the graphical corrections we have calculated, we obtain 
the following RG flow equations: 



dB{e) 
di 



= d-l 



2X-Y633 



(136) 



^^(d-5 + co + 2x+^g3]K, (137) 



M =H-3 + - + 2x+^53 



de 

dD(£) 
d£ 



d-3 + 3x- 



32 



53 



32 



54 g 



(138) 



(139) 



d-3 + uj + 2x+-^g3- ^94] D 



+ T^^ (94 - 53) , 



dwji) 
di 



16 



= d-5 



(140) 



3.9i 



3254 



3 15 

g53-^54)^. 



(141) 



These equations are not very transparent telling the 
physics since they depend on the rescaling factor uj and 
X, which are artificial. Instead, we should focus on the 
flow of the dimensionless couplings. Combining the above 
flow equations with the definition of 53 and 174, we find 



^33 . 9 

-dF = ('"i6-^^ 



dg4 
di 



15 
32 



94 



13 

-^93)93 : 

13 



32 



33 



£5i 
32 34 



)g4 



(142) 
(143) 
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where e = 5 — d. 

Obviously the Gaussian fixed point = 0, g| = 0) is 
one of the fixed points of Eqs. p42p and (|143p . To check 
it's stabihty, we hnearize the two equations around it by 
writing g3 — + Sg^, 34 = + Sg^, and get, to hnear 
order in (5(73,4, 



reduces to 



M 
dSgi 

d£ 



<^^g3 : 

eSgi. 



(144) 
(145) 



These two hnearized differential equations tell us that for 
e > {d < 5), given a small departure from the Gaussian 
fixed point, both (73 and (74 will leave it exponentially, 
which implies that the Gaussian fixed point is unstable 
in the physical dimension d = 3. 

In addition to the Gaussian fixed point, these RG ffow 
equations (|142[ 1143)) have other two fixed points. One of 
them is exactly the one found in the isotropic disordered 
smectic A problem [^], for which the Hamiltonian is ro- 
tation invariant. This is not a coincidence, because the 
Hamiltonian we start with is more general in the sense 
that it has no symmetry bound (g/B)^ = w/B, so our so- 
lutions must include this specific fixed point as a special 
case. To find this particular fixed point , we can en- 
force this symmetry condition on the RG fiow equations 
(|1421 I143P by setting 33 — g^. After this enforcement 
both of those two equations reduce to the recursion rela- 
tion found by Radzihovsky and Toner Q for the isotropic 
problem: 



dgs 

de 



which leads to 53 = 32e/5. However, this fixed point is 
not stable in our more general problem. We can examine 
this instability by using a similar linearization analysis as 
for the Gaussian fixed point. Writing (73 = 32e/5 -t- Sgs, 
(74 ~ 32e/5-f (5(74 and plugging them into equations (jl42[) 
and (|143p . we obtain 

dSgs 13e 18e 

^ = — ^53-— <5ff4, (146) 

d6g4 7e 12e 

^ = jS9s-—Sg„ (147) 

which has the solution near the vicinity of the fixed point: 

Sgs = Cie-'^ -I- 18C2e6'^/^ 
,554 = Cie-'^ -I- TCse^^^/^ 

The Ci and C2 in the solution are constants determined 
by the initial condition ((5,g3)o and (694)0. This solution 
shows that in general, this fixed point is also unstable, 
although it is stable given that {Sg3)o = {694)0, which is 
just a verification of its stability in the isotropic disor- 
dered smectic A problem . 

The third fixed point can be obtained by assuming that 
(73 fiows to zero. Under this assumption equation (jl43p 



dg4 , 15 

^ = - 5^.94)94 , 



(148) 



which has a non-Gaussian fixed point g\ = 32e/15. The 
linearized differential equations around this fixed point 
are 



di 
d6g4 



1 



13 

df ~ 15 
which have the solution: 



:<^6g3 , 
e6g3 - e6g4 



(553 = Cle-^^/^ 
13 
12' 



^54 = ^<3ie~"^/^ + C2e" 



(149) 
(150) 

(151) 
(152) 



Note that both eigenvalues are negative for d < 5, which 
implies that this fixed point is stable and hence controls 
the second order phase transition. 



Ratd-Toner Fix Point 




FIG. 20: Graphical RG flow of the dimensionless coupling 
and (?4. Point a is the stable fixed point (0, 32e/15) 



IV. CRITICAL EXPONENTS 

The RG we derived in the previous section is only valid 
in the critical region, and hence named "critical RG" . If 
the system is right at the critical point, the critical RG 
holds for arbitrary £, otherwise it breaks down at i* which 
is determined by D{£*) = K{£*)A'^. In this section we 
first identify the critical point, then study how fast the 
critical RG breaks down if the system is slightly off the 
critical point, so as to calculate the critical exponents. 
For convenience, we choose a special rescaling 



uj + 2x = 5-d- 



1 

32 



53, 



(153) 
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which fixes K at its initial value Kq. The flow Eq. (|140p 
then becomes 



d£ 



2 + \9,-y.]D 



+ -r7:Ko{gi- gs) 
16 



(154) 



Right at the critical point D{£) also flows to a flxed value 
D*, which can be calculated by plugging = 0, gl = 
32e/15 into the Eq. (fTM)) and requiring dD/d£ = 0. 
Thus the fixed point controlling the phase transition is 
identified in I?-(73-(74-space as 



= 0, 54 = 32e/15. 
D* = -e/5. 



(155) 
(156) 



For initial values Dq, g^ and 174 off from the critical plane, 
D{i) departures exponentially under the RG and for large 
£, behaves as 



D{e) oc C3e^"' cx (T - TAc)e^'' 



(165) 



The system is in the A phase for T > Tag (i-e., C3 > 0) 
and in the C phase for T < Tag (i-e., C3 < 0). The 
universality classes of both phases will be discussed in 
section I VII 

Now we calculate the critical exponents , Vz- These 
two exponents are defined via the two correlation lengths 
and £,z by equations 



6 cx \T-TAcr'- 



(166) 
(167) 



Linearizing the fiow Eqs. (|142p . (|143p and (|154p around 
this fixed point, we obtain (I149p . (I150p and 



dSD{£) 
di 



Xd6D +{\d*- ^Ko ] 6g3 



where 



5 



(157) 



(158) 



The solutions of these flow equations are given by (|15ip . 
(HHl) and 



6D{i) = C3e^" 



128 ' 



ee/5 



iKoC2 

32 



e-^^ (159) 



where Ci_2,3 are variables determined by the initial con- 
dition via 



Cl 
C2 



.93" 

9l 



53. 

ff4 



C3 = Do-D* 



128 



53 J 



32 



(160) 
(161) 

(162) 



where Dq, g^ and g^ are the "complete bare" values of D, 
gz and 54, respectively. According to the solution given 
by (|159p . the flxed point is unstable with respect to D 
unless C3 = 0, which deflnes exactly the critical point. 
After some algebra, the critical point is identifled as 



AC) 



32 ^3 



32 



(163) 



which actually describes a plane in D-g^-g^ space. This 
shows that the true critical point has dependence on ^3 
and gl, and is not exactly at Dq{T) — 0, which is pre- 
dicted by the mean-fleld theory. C3 measures the devia- 
tion from the critical point and is often assumed to be 



C3 cx T - Tag ■ 



(164) 



These two correlation lengths are the characteristic 
lengths beyond which the critical behavior crosses over to 
either the high temperature (smectic A) or low tempera- 
ture (smectic C) behavior. More speciflcally, the correla- 
tion length satisfles K{£^±)^'^ = D, which, under the 
RG, is mapped to 



with 



(168) 



(169) 



where due to the special rescaling (I153p . K{i*) is fixed 
at Kq and D{£) is given by (|165p . Thus we obtain 



r cx--^in|r-r^cl, 

Ad 



(170) 



which combined with Eq. (|169p . gives 

a«|r-rAcr/^". (171) 

Therefore, the RG predicts 

Re 

(172) 



iy± = - = - -\ hO (e^ ) . 

A 2 20 ^ 



The other correlation length can be obtained using the 
scaling law cx this leads to 

3e 



1 



10 



(173) 



In order to calculate the specific heat exponent a, we 
need to find the temperature-dependence of the free en- 
ergy. Since the free energy F is invariant under the RG, 
we can write it as 

F = Vrfr{Kit),B{e*),Din,A{n,9{n,^{n) 

where Vr , fr are the renormalized volume and free energy 
density, Vq is the physical volume before renormalization. 
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To factor out the temperature dependence of fr as much 
as possible, we make a special choice of £*, oj and x- 



t = In(aA), 



(174) 
(175) 

(176) 



This choice fixes K and B so that they are ^-independent 
and hence temperature-independent. The choice of £* 
makes D{1*) = KqA?' and hence also temperature inde- 
pendent, and for such large D{£*) fluctuations are small, 
and fr can be evaluated by ignoring the anharmonic 
terms. In spite of all the advantages just described, 
fr is a function of A{£*), which still has temperature- 
dependence. Fortunately by taking dfr/dA{£*), which 
can be easily calculated, we find fr is a linear function of 
A(^*). Finally the free energy can be nicely written as 

F = Voe-'^^-^'^^'-^^' {A{£*) + constant) 



oc \T-Tc 



X 



IT - r,r^(-'+''^-'") + constant 



(177) 



Now taking d^F/dT^ and keeping the most divergent 
part, we get the specific heat exponent 



1 H Vt 



(178) 



V. ANOMALOUS ELASTICITY RIGHT AT THE 
CRITICAL POINT 

In this section we show that the existence of the stable 
fixed point leads to anomalous elasticity by using tra- 
jectory integral matching. For the sake of simplicity we 
assume that the system is right at the critical point, so 
that we do not need to worry about the crossover from the 
critical region to the A oi C phase, which will be taken 
care of when we discuss correlations functions in section 
IVIIII Let us cal culate the d isorder averag ed correlation 
function G{q) = {\u{q)\^) - {u{-q)) {u{q)). The RG es- 
tablishes a connection between a correlation function at a 
small wavevector (which is impossible to calculate in per- 
turbation theory due to the infrared divergence) and the 
same correlation function at a large wavevector, which 
can be easily calculated in a controlled perturbation the- 
ory. This connection for G{q) is 

G(g, B, K, At,g, w) ^ e(rf-i+-+2xK x 
Giq^c', q,c^\ B{£),K{£), At{£),g{£),w{£)) , 

(179) 

where the prefactor on the right-hand side comes from 
the dimensional and field rescaling. First let us consider 



the special case Qz = 0. The rescaling variable £* is cho- 
sen as 



q±e 



A. 



(180) 



We also choose q± sufficiently small such that 53 (^*) and 
g4{£*) has reached the stable fixed point = 0, gl = 
32e/15. Eliminating £* in favor of q±, we get 



Giq^,qz = 0,B,K,At,g,w) = 



A 



{d-l+u+2x) 



xG(A, 0, B{£*),K{£*), At{£*),g{£*), w{£*)). 

Evaluating the right-hand side using a perturbation the- 
ory, we obtain 

G(q^,qz = 0,B,K,Aug,w) 



1 



K'^K{£*) 
1 



A 



{d-l+uj+2x) 



K{quqz=0)qi ' 

where 

K{£*) = K,{^hL^f-''^''^^''{^hLqi 



(181) 



-(d-5-l-w+2x-S3/32) 



which is obtained by integrating the flow Eq. (|137p . 
and the anomalous bend modulus which diverges at long 
length scales is defined as 



K{q^,qz=i)) = K^{q^^j,L)-^'^, 
with an anomalous exponent 



(182) 



(183) 



^]^^ is the nonlinear crossover length along _L directions, 
beyond which the perturbation theory breaks down. 

The above calculations can be generalized for arbitrary 
direction of q by using a more sophisticated choice of £* , 
which satisfies 

K{t)K^ = K{t){q^e''f + B(r)(g,e"^*)2 . (184) 



This choice ensures that the right-hand side of (|179p can 
be evaluated in a controlled perturbation theory and is 
consistent with (|180p in the special case qz — 0. Solving 
P84p we obtain 



NL 



where the anisotropy exponent C, is given by 



(185) 



(186) 



rjB is defined as the anomalous exponent of the compres- 
sion modulus B 



5(9^,9. = 0) = So('7±ewi)"", 



(187) 
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and given by 



VB = Y^53 



fe{x) is a scaling function 



1, X < 1 

X^/'^, X > 1 



(188) 



(189) 



Plugging (|185p into p79p leads to the general q- 
dependence of both K and B, which is summarized in 

The g-dependence of At can not be cal c ulated by using 
the correlation function G{q) = {\u{q)\'^) — {u{—q)) {u{q)) 
since it is independent of At. Instead we perform tra- 
jectory integral matching on the correlation function 
{u{q)){u{—q)) = Atq^G{q)^ . After a similar calculation 
for the special case = 0, we find At also diverges at 
small q± as 

At(qx, = 0) = AO(qxC]i^L)-'" , (190) 

where the anomalous exponent rjt is given by 

Vt = ^33- (191) 

Generalizing the calculation for arbitrary q, we get the 
result summarized in p2p . 

Since 53 is zero, rjK.t sue zero to 0(e), and rjB is zero 
exactly. To calculate rjj^ t to O(e^), we performed a two- 
loop RG calculation showing that 



^t CAe' + 0(e'), 



(192) 



where Cr = (32 In (4/3) - 10)/225, Ca = (12 In (4/3) 
l/3)/225. 



VI. A AND C PHASES 

In this section we will identify the universality classes 
of both the high temperature (A) and low temperature 
(C) phases. For now let us assume that the system is 
far away from the critical point so that we do not need 
to worry about the crossovers between the critical re- 
gion and the two phases, which will be discussed in sec- 
tion IVIIII First we discuss the A phase. The model 
we start with is the Hamiltonian ((80)l . Since at long 
wavelength the tilt term Do\V ±u\'^ /2 dominates the bend 
term K{W'^u)'^ /2, the latter can be neglected. The model 
reduces to 



H = d 



V I u 



(193) 



For this model a simple power counting shows that both 
the anharmonic terms and the random tilt disorder are 
irrelevant in d > 2. Thus in d = 3 the effective model for 
the A phase is simply 



H 



B 



D 
"2" 



V I It 



Vp{u-c^{r)) , 
(194) 



which is just the random field XY model with anisotropic 
stiffness. This model has been studied extensively, and 
the correlation {{u{r) — u(0))^) has a logarithmic diver- 
gence as r — > 00. As a result the density correlation 
function for the A phase decays as a power law, which 
implies that the A phase only has quasi-long-range trans- 
lational order. However, unlike the conventional smectic 
A phase, in which the power law exponent depends on 
the temperature, here the exponent is universal. These 
correlation functions will be calculated in detail in section 

Enn 

Deep in the C phase, Z3|V_LUp/2 is negative. Let us 
first use a mean field theory to find the new disorder 
averaged gvonnA state of the rfisordererf Hamiltonian (|80p . 
Since the two disordered interactions in the Hamiltonian 
are totally random and isotropic, they should not affect 
the mean field theory. Therefore we ignore them for the 
moment and work on the clean Hamiltonian 



H = d'^ 



V I u 





4" 







(195) 



which can be rewritten as 
'B 
"2 



H 



d'^r 



^1 
2B' 



K 



2D 



w' 



2w' 



where 



,' = w^ (g^/B) . 



(196) 



(197) 



w' = is the tricritical point separating the first order 
{w' < 0) and second order {w' > 0) phase transitions. 
Assuming w' > since we are interested in the second 
order phase transition, the ground state of the Hamilto- 
nian (jl96p is given by 



2D 



Bw' 



(198) 
(199) 



Now we go back to the disordered Hamiltonian ([50)1 and 
expand it around its disorder averaged minimum by mak- 
ing the substitution 



2D 



^ 9D ^ , 
Bw' 



(200) 
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where we have assumed that in the C phase the sys- 
tem spontaneously breaks azimuthal rotation symmetry 
and the layer normal tilts along x, which is an arbitrary 
direction within _L plane. Also for simplicity we throw 
away the random field term by using a posterior reason- 



ing; that is, if the C phase is "m = 1 Bragg glass" , then 
the random field disorder is irrelevant [8|. For brevity, 
we define 6*0 = \J—2D/w'. After the substitution and 
simplification, the Hamiltonian ((80)) becomes 



H\u'] = ^ I d'^r 



+ wdo{dxu') 



w 



h (r) • Vu' 



(201) 



where s denotes the direction which is perpendicular to replica trick, we get 
both X and z. Treating the tilt disorder by using the 



H\u' 



K (V^O' + B {d^u'J + wOl {d^u'S - 2(?0o(5.O(9.<) - g{d,u'^) 



-I- w6'o(9j;U^) |VsUq| + — 



2T 



(202) 



a,/3=l 



Q,/3=l 



whose universality class is still not clear. Then we make with 
the following coordinate transformation: 



rz' = {g/wdo)r^ + r^, 



(203) 



which, in momentum space, corresponds to the transfor- 
mation 



w V —2D 



(205) 



qz' = Qz , 

Qx' = Qx- Tqz , 



(204) 



After this coordinate transformation the Hamiltonian 
pl^ becomes 



^K] - o h'^-'T. K{wi,u':)^+B{dz:u'j+^{d,.u'j+m<) 

•' a=l - 

, „ n 

- ^ / E [As' {y.<-ys'u'p) + A,,{dW^){d,,u'^,) + Az,{dzK){dz,u'^) 

a, 13=1 
If " 



(206) 



Q,/3=l 
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with 



B 




(207) 




w 


k 


= K, 


(208) 


7 




(209) 


9 


= w9q , 


(210) 


As' 


= A., 

' — " ) 


(211) 




= At, 


(212) 




2gAt 


(213) 
















(214) 



It is easy to check that 7, g and w satisfy the magic rela- 
tion g = Y^wFy. This is not a coincidence but due to the 
symmetry, that is, the Hamiltonian must be invariant un- 
der the rotation about z-axis since the environment (i.e., 
the aerogel stretched along z) is azimuthally symmetric. 
Because of this symmetry the model (|206p belongs to the 
universality class of "m — 1 smectic Bragg glass" . 



VII. WAVEVECTOR-DEPENDENCES OF 
IRRELEVANT DISORDER VARIANCES 



In our RG calculations in section IIIII we have not in- 
cluded in the Hamiltonian ([55)1 the random compression 
disorder 



Ac /■ " 



(215) 



A simple power counting shows that this term is irrele- 
vant and thus has no effect on the critical behavior and 
anomalous elasticity. However, while it does not affect 
the anomalous elasticity, the disorder invariance Ac itself 
develops strong power-law dependence on the wavevec- 
tor at long wavelength. Now we calculate the power-law 
exponent, which we denote as rjc, to 0(e). This is nec- 
essary since the random compression disorder has non- 
trivial contribution to the correlation function (|u((7)p)j 
and therefore affects the light scattering predictions. 

The most important one-loop graphical correction to 
Ac comes from the Feynman diagram in Fig. [71 from 
which we obtain 



^Ac 



2tt 



i2n) 



d-l 



5 

64 



9 



{Kql+Bqlt 
A'^-\ (216) 



Thus using the dimension and field rescaling we used in 
section IhII the RG flow equation of Ac to one-loop order 



is given by 

rfAc(^) 



1-UJ + 2X+ —35 A 



(217) 



de \ 64- 

where the dimensionless coupling g^ is defined as 

55 = BAtg^/KA, . (218) 
Combining this RG flow equation with the flow Eqs. 

inm \m\ \mi nmim . we get 



dg5 
M 



2- -e 
5 



64- 



55 35, 



(219) 



which flows to a stable nontrivial flxed point g^ — 
64 (10 — e) /25. Then the wavevector-dependence of 
Ac can be calculated by performing trajectory integral 
matching on the correlation function {u{—q)){u{q)) = 
{Atq'^ + Acql)G{q}'^ . The calculations are essentially the 
same as those in section IVl and will not be repeated here 
again. The result is given in ([5^ with the anomalous 
exponent 



Vc 



5 , „ 1 
755 = 2 - -e 



(220) 



64-^ 

Note that rjc is nonzero even to the zeroth order of e, 
which is quite common for irrelevant variables. 

It is also interesting to point out that rj^ is not fully in- 
dependent, but are related to other anomalous exponents 
by an exact scaling relation, which is implied by the fact 
that (75 flows to a non-zero stable fixed point. For large 
enough £, 35 reaches the fixed point and thus 

din 35 



M 



= 0. 



(221) 



This equation, after decomposing 35 into small pieces by 
its definition (|218p . leads to 

d\nB dhiAf dlnga dlnA' dhiAc 



dl 



d£ 



de 



dl 



de 



0. 



(222) 



After plugging the flow Eqs. ([131 [T371 [HHl [B3l (HIl) into 
the above equation, the rescaling factors x ^-nd uj vanish 
and we are left with 



2 + f]t - Vb - Vk - V3 ■ 



(223) 
rjx and 773 to 



Since we already know the value of rjt, Vb, 
0(e), using ([^ also leads to ([^ . 

Likewise in the model of the C phase, which is given by 
(j206p . there are also disordered terms which are irrelevant 
in the RG calculation but important for making light 
scattering predictions. These terms are 



A^ 
2T 

2T 



d'^r 



E 

a, (3=1 



rY.dz 

a, 13=1 



A 

It 



^ d' 



a, 13=1 



•5, 
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Az' and ^x'z' have no wavevector-dependences since 
there are no graphical corrections to them. 

As for Aa;', we wiU use a different RG, namely the 
RG calculation for "m = 1 Brag glass", to show the 
power-law wavevector-dependence of A^^/ and calculate 
the power-law exponent, which we denote as f]^' ■ For 
convenience, we transform the Hamiltonian (j206p into 



the exact form studied in [^. This transformation is ac- 
complished by doing the field rescaling 

v! = Qqv!' . (224) 

In terms of u" the Hamiltonian (|206p can be written as 



1 /■ " 



K' (V2,<)' + B' {d.,ulf + y {d.<f + 7'(5.'<) 



2 ryl 



-. „ n 

J E [A;'(v.'<-v,m;^)-HA;,(a.,<)(a.m^) 



Q,/3 = l 



(225) 



where the coefficients of the quadratic terms are related 
to the original ones by 



s' ,x' ) 



(226) 



For this model, the detailed RG calculation can be found 
in Q. Here we describe the calculation very briefly. A 
hyper-cylindrical Brillouin zone is used: —A < g^/ < A, 
— oo < < oo, < \qz' I < oo, where A is an ultra-violet 
cutoff. Note that based on the arguments in the previ- 
ous section, A has to be replaced by when the sys- 
tem is within the critical region. The displacement field 
is separated into high and low wavevector components, 
= ""a"^ (r) + ■"q^ where u'^^ (r) has support 
in the hyper-cylindrical shell Ae~^ < \qs'\ < A, — oo < 
Qx' < oo, < \qz' \ < oo, and (r) has support in the 
remainder of the hyper-cylinder (i.e., < Ig^'l < Ae^^, 
— oo < qx> < oo, < \qz'\ < oo). Then integrate out 
the high wavevector part u^^ (r), and rescale the length 
and long wavelength part of the fields with 



to restore the UV cutoff back to A. 



u'^{r') so as 



The underlying symmetry of the model (j225p guaran- 
tees that the unit {dx'u'^ — (Vs'U^)^/2) is graphically 
renormalized as a whole. Therefore it is convenient to 
choose the rescaling which also preserve the unit; the ap- 
propriate choice is x' = 2 — ujx' ■ 

Using the hard continuation (namely, keeping the num- 
ber of the soft directions fixed at 1 for all spacial dimen- 
sion d) , the RG to one-loop order gives the following flow 



equations: 

M 

dK'je) 
d£ 



5 -I- (d — 2)ujzi — ?>ijJxi - 
l + {d- 2)ujz' - Ux' + 

3 + (rf — i)uJz' ~ ^x' + 



356 



32\/2 



.96 



8\/2. 

56 



32V2. 



7 , 
(227) 

K' (228) 
(229) 



where 56 = A'^.{i / {K''^"^ B"^-'^)f^Cd-ih-^'^-'^ is a di- 
mensionless measure of disorder, whose flow equation is 



2e56 



15 



:56^ 



(230) 



where e = 7/2 — d. This equation has a stable fixed point 



56 



128^2 



(231) 



Using trajectory integral matching we can calculate the 
wavevector-dependences of 7', K' and A^, , which, com- 
bined with (|226p , leads to the wavevector-dependences of 
K, As/ and 7 summarized by ([H [71 H]), with the expo- 
nents 



Vk 

Vs' 

1% 



= 2- 



+ Vk 



9_ 

2 ' 



56 



8V2 15 
32V2 15 ^ 



(232) 
(233) 
(234) 

(235) 
(236) 
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That both ge and gy flows to non-zero fixed points im- 
pUes the exact scahng relations ^T8\ fT9l ). For the deriva- 
tion, please refer to how we derive (|223p . 



In addition, using soft continuation (i.e.; keeping the 
number of hard directions fixed at 2) obtains 



2 + i + 0{f) 



(241) 



and the same exact scaling relation given by (|19p . where 
e = 4 — d. The numerical estimate and error bar for 
r]x' can be obtained using the approach described in Ap- 
pendix El 



VIII. CORRELATION FUNCTIONS 



FIG. 21: The graphical correction to A'^,. 



So far we have briefiy reviewed the RG calculation 
for "m — 1 Brag glass", now we use it to calculate the 
wavevector-dependence of the irrelevant disorder invari- 
ance A^, . The Feynman diagram for the one-loop graph- 
ical correction to Af^, is presented in Fig. [2TJ An analysis 
of the diagram gives 



SA' = 



2ti 



dqs 



„8 



21 



128\/2 



(2^) 



.97 : 



B'qi, 



I'll 



(237) 



where 177 is another dimensionless coupling defined as 
57 = g67'A;,A-V(iv:'A^,)- Thus the RG fiow of A^, 
to one-loop order is given by 



dl 



5 + {d ~ 2)lUz' — SuJx' + 



21 



128V2 



97 



(238) 



Combining this flow equation with ()227p . (I228|) and 
(|229p . we obtain the flow equation of the dimensionless 
couphng 57: 



dgi 
M 



8 



21 



57 



128%/2 



97 



(239) 



which flows to the fixed point = 256A/2(5-4e)/105. At 
this point, our experience from the previous calculations 
immediately gives us the wavevector-dependence of A^,, 
which, combined with (|226p . leads to the wavevector- 
dependence of Ax' summarized by The power-law 
exponent 77^:' is given by 



21 



Tlx' 



128%/2 



57 



(240) 



In this section we discuss correlation functions and 
make X-ray and light scattering predictions. We start 
with the calculation of the disorder averaged fluctuations 
in momentum space <\u{q)\'^ >. In our problem this 
quantity has three main distinct types of behavior. They 
are characterized by the fluctuations in the A phase, the 
fluctuations in the C phase and the fluctuations right at 
the critical point, respectively. The first two types have 
been well studied and known, and the third type is unique 
to this problem. Our main interest is the rich crossovers 
between these distinct behaviors near the critical point, 
which leads to amazing experimental consequences. We 
warn the reader of heavy algebra in this section. 

As we pointed out earlier, in some cases while certain 
type of disorder is irrelevant in the RG sense, it can have 
important contributions to correlation functions. Thus 
we prefer to use the complete model given by (|82p to 



calculate < |w(q)P >. 

For large qa (i.e., q±_ <C CJ^ or q^ <C ^^^), the model 
can be treated by the critical RG, and we get 



< \u{q)\^ > = < \u{q)\^ >AC 

At (q) ql 



{Bql + D{T,q)ql+K{ci)qiY 

^c{q)ql 

{Bql + D{T,q)ql+K{q)qif 



(242) 



where K and Aj.c are all g-dependent and given by (j271 
SH), B is g- independent, D has dependences on both q 
and the temperature: 



D{q,T^TAc) 

(T - TAcUhL^l-T"^ a Liz « (e^L9±)'^ 

(243) 



where 77/5 = 2—rix — l/iy±- The first and the second terms 
in the Eq. (|242p come from the random tilt and random 



31 



compression disorders, respectively. In this case the ran- 
dom field disorder is not only irrelevant in the RG calcula- 
tion, but also has negligible contributions to < > 
in any region in q'-space, and is thus neglected. 

For small qs (i.e., q± <C or <^ Cj^), the tilt term 
in the model dominates the bend term, and the critical 
RG is no longer valid. A different RG needs to be used to 



treat the model. First let us derive the effective model for 
small qs. The derivation has two steps. The first step is to 
run the critical RG to the crossover point £* = In(A^x) to 
obtain a model for the rescaled system. The second step 
is to undo the dimension and field rescaling we performed 
during the RG. implementing this procedure we obtain 



K{T) (Vlua) -1- B {d.u^Y ~ g{T){d,Uo.) (V^u^) 



D{T)\V^u^\ 



E 



V_LMa • V_lW^ H ^l^O^Ua ■ d^Up + —Ap{Ua - Up) 



(244) 



where the coefficients are renormalized by the critical 
fluctuations and hence temperature-dependent. 

In the A-side critical region, based on the argument 
we give in section IVIl the model (|244p belongs to the 
universality class of "random field XY model" . We treat 
it using the functional RG 3 and obtain 



< |u(g)|2 > = < \u{q)\^ >A 



At {q,T)ql 



{Bql + D{T)ql+K{T)q\_f 



^c{q.T)ql 



{Bql + D{T)ql+K{T)q\:) 
CB^D 

{Bql + D{T)ql+K{T)q\)^ql 

(245) 

where the third piece comes from the random field disor- 
der. The random tilt and compression disorder are irrele- 
vant in the RG calculation; however, their contributions 
to < |u((7)p > are important in some region in q'-space 
and hence kept. The constants D, At.c in (|245p are 
not wavevector-dependent and given by 



B 
K 
D 

Atr 



Boi^i^jUr- ^Bo; 



In the C-side critical region, we have to expand the 
Hamiltonian (|244p around the disorder-averaged mini- 
mum. Following the procedure in section IVIl leads to 
the model (|206p . where for convenience a special coordi- 
nate system other than the lab one has been used. The 
relations between these two coordinate systems are con- 
trolled by the parameter F and, in Fourier space, given 
in Eqs. ([51 [51 flU)) . F is now renormalized by critical 



fluctuations and hence temperature-dependent, given by 

'iK+'ia 



(aA)'- 



(246) 



The exponent is equal to 4/5 in d = 3 to 0(e); there- 
fore, F is expected to be very large as T ^ Tag- "^^^ 
derivation of this result is given in appendix [Bl The 
temperature-dependences of the coefficients in the model 
P06p can be derived in a similar way and are given by 



7. 



Kc = KoiC^J^^r'- oc (Tag - T)-"--- ,(247) 
Do{^nl/UY'' « {Tag - r)^^-^).. ^ 

(248) 

A?(^^L/e±)-'" « {Tag - T)"""-- , (249) 



A 

<x {T-TAGr"^ 



6 



(250) 



These coefficients have been named "half-dressed" values 
(i.e. unrenormalized by the C fiuctuations) in Eqs. ^ 
[71 [6]). Treating the model (|206p using the RG for "m = 1 
smectic Bragg glass" Sj , we obtain 



< \u(q)\^ > = < \u{q)\^ >c 



As' {q',T)ql 



^{q',T)ql,+Bql+k{q',T)ql, 

A.^ (g^^)e 

(-l{q',T)ql,+Bql,+k{q',T)q^^, 
A.' {T)ql 



l{q'.T)ql,+Bql,+K{q',T)q%^ 

(251) 

where the g' -dependence's of if, 7, As'^x' are given in 
Eqs. ([ni[71[S]), B and A^/ are -independent. The first 
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and the second terms correspond to the contributions 
from the random tilt disorder along the soft and the hard 
directions in the original "m = 1 smectic" problem, re- 
spectively, and the third term corresponds to the contri- 
bution from the random compression disorder in the same 
problem. The random field disorder is neglected since its 



contribution to < |u((7)p > is always subdominant to the 
others. 

Thus, near the critical point the behavior of 
< > can be summarized in the following compact 

way: 



r < \u{q)\'^ >A, g± « ei\ « for T> Tag 

< \uiq)\^ > = I < \u{q)\^ >AC , q± » Cj' or » C\ ■ (252) 

[ <\u{q)\^>C, q± « qz « for T < Tac 



The X-ray Scattering pattern for both A and C phases 
has been studied before. It is quasi sharp, isotropic for 
the A phase, broad, anisotropic for the C phase. In the 
following we will discuss the crossover between the two 
distinct patterns in the critical region. The X-ray scat- 
tering intensity is given by the Fourier transform of the 
thermal and disorder averaged p — p correlation function: 

/(g) oc J d'^r'd'^r" (p(f ')p(f "))e-*«-(""-"") , (253) 

where p{r) is the molecule density. p{r) can be expanded 



in a Fourier series with period a, the distance between 
nearest layers, via 

oo 

p{i^= p„e"9o(z+«(r)) ^ (254) 



where p„ is the (complex) amplitude of the n'th harmonic 
of the smectic density wave, u{r) is layer displacement at 
the point f. Inserting the decomposition (|254p into (|253p 
gives 



/(g) cx 



E 



PnPn 



d^'d^r" 



^-iq-(r'-r")^iqa(i 



(255) 



Changing variables of integration from r ' to the differ- 
ent variable r = r ' — r" and using the fact that, in a 
homogeneous system 



with 



{^giqa{nu{r+r")+n'u{r"))^ — (^g^iqa(nu{r)+n' u{a))^ (256) 

the integral over r" now just gives V5n+n', where V is 
the volume of the system and (5„+„' is a Kronecker delta 
which can be used to collapse the sum on n' to the single 
term n' — —n. Doing so, we obtain 



I{q) cx Vj2\pn? / d^'r e'("9"^-9^-^>„(r") , (257) 
where 

F„(f) = (^e**(™(*')+"'"(o))^ . (258) 

Now approximating the u - fluctuation as Gaussian, we 
can write 



^^„(f)=cxp(-^C(f) 



(259) 



C(f) ^ ([u{r)-u{0)f) 



d'^k 



2[l~cos{r-q)]{u{q)u{-q)). 



(260) 



Putting l[259l) into ([2571) we get 



Qo p / 2 2 

/ (g) (X ^ d'^r e*("'?°^-«^-^^exp ( -!i-^C(f) 



(261) 

For q near the Bragg peaks at the mqoZ , the sum is read- 
ily seen to be dominated by the n = m term; thus, the 
scattering near the Bragg peaks is essentially the Fourier 
transform of exp(— n^goC(r)/2). 

Right at the critical point (T = Tac), in spite of the 
symmetry broken we mentioned earlier, the fluctuations 
of the system is still qualitatively like those of "Glassy 
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smectic A", hence the X-ray scattering pattern is also 
broad and anisotropic. 



For T TJc, plugging Eq. (|252)) into Eq. ((260l) and 
performing an asymptotic analysis, we obtain, in 3D, 



J 



«6 



<6 



r± > or > 



where F = 2 + — (3?7i<-/2), fr is another universal 
scaling function, « 1.10, g = (X/UM^/^nl)'"'^^- For 
small r, i.e., r± ^ C-Li ^2 ^ Cz, C(^) increases as a power 
law, while for big r, i.e., r± ^ or rz ^ ^z, C{'r) in- 
creases logarithmically. This indicates that Fn{r) varies 
faster (exponentially) for small r {r±_ <^ ^±,rz ^ ^z), 
slower (as a power law) for big r {rj_ ^ or 3> 
^z). So the X-ray scattering intensity I{Sq) for large 
= \G — q\, which is dominated by {5q), defined 
as the contribution coming from the fast varying part 
of Fn{r), should be broad and anisotropic, qualitatively 
like a Lorentzian squared. The line widths Sq^, 6q^ of 
this Lorentzian squared are defined as Sq^ = (^f )~^ and 
Sq^ ~ (CI) "'^ with and satisfying 



C(0,Cz") = a2, 
C(ei,0) = a2. 

Solving these two equations we get 

2/r 
2C/r 



(262) 
(263) 



- (- 



(264) 
(265) 



The temperature dependence of and could be used 
to determine the exponents 2/F, 2C,/ F, rjk, and 77A since 
the bulk K{F) and B{F) that implicitly appear in Eqs. 
(|2641 1265P have temperature dependence that can be ex- 
tracted from measurements on bulk materials. 

For small 6q, I{5q) is dominated by I^{5q), which is 
defined as the contribution coming from the slow varying 
part of Fn{r). Thus, I{5q) is quasi-sharp and isotropic, 
and diverges as a power law as 5q approaches 0. 

The crossover between these two distinct scattering 
patterns is defined as the point where the two contri- 
butions become comparable, that is, P{5q) — P{5q). 
Assuming that this happens at very small \5q\, so that 
I^{Sq) can be treated as a constant {^^) . We will ver- 
ify this assumption as a posteria. Let us first calculate 



the crossover in _L directions 

I{Sq^ = Sql,Sqz^O) OC (5q^)-3+0.55n= ^-0.55n^ 



xe 



which gives 



Sq'l cx |r - Fac 
X exp( 



|0.55ri^O/(3-0.55n^) 



(266) 



iF-FAcr""). (267) 



3 - 0.55n2 



where Q — Fv±. 
found to be 



Likewise the crossover in z direction is 



\T - Fac 



(268) 



Both 8q\ and iq'i, become extremely small when F 
Ftr^, which is consistent with our earlier assumption. 



As for F —^ F^Q, since the X-ray scattering pattern is 
always broad, there is no significant crossover which can 
be detected by experiments. 

Now we investigate light scattering behavior near the 
critical region. The light scattering intensity is propor- 
tional to a linea r sum of the fluct uations of the nematic 
director Aij{5n^{—q)6njr{q)) [l2|, where the values 
of Aij are determined by the electric polarization direc- 
tions of the incident and transmitted light. Since in Sec. 
|TT]we have shown that the fluctuations of fi and the layer 
normal N are bound together, the light scattering inten- 
sity is thus proportional to a linear sum of Cij{q), which 
is defined as 



a,iq) ^ {N,^{-q)N^{q)) = L^(-Z>i< Hq)\' > . 

(269) 

We now derive the eight distinct regions with qualita- 
tively different wavevector-dependences of Cij{q) in the 
A-side critical region, which are illustrated in Fig. 01 For 
big g's (i.e., q±_ ^ or qz » which corresponds to 
the regions outside the rectangle O JFI ) , Cij (q) is given 
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by 



{Bql + D{T)ql 
Ac {q)ql 



K{q)q\f 



{Bql + D{T)ql + K{q)qiy 



(270) 



where in the denominator Dq\ is neghgible compared 
to Kq\. The fmiction for the crossover (i.e., locus EG) 
between the random compression and tilt can be obtained 
by doing 

At {q)qi 



{D{T)ql + Bql + K{q,T)qi) 

Ac (g) qlq\ 

{D{T)ql + Bql + K{q,T)qi) 



(271) 



which leads to 



(272) 

Using the scaling relation (|223p . this equation can be 
further transformed into 



qzi,_ 



NL 



q±i 



NL 



A?/AO 



(273) 



where A = C/(l +'73/2). Above EG (i.e., in region 8), 
Cij (q) is dominated by the random compression fluctua- 
tions and given by 





fA,qlql\ 
V B^qt ) 


4(^1 1 




L^q) 1 


( A^A" (q^ 





7^ - {q^CkL) 



(274) 

where we remind reader that A is the smectic penetration 
length defined by A = Kq / Bq , and ^^'^ are the non- 
linear crossover lengths defined in Eqs. jSOlISl]). Below 
EG, Cij{q) is dominated by the random tilt fluctuations 
but still has two different wavevector-dependences due 
to the crossover (locus FH) between Kq\ and Bq^ in the 
common denominator, which satisfies 

q^^NL = {q^^NLf (275) 

and is below locus EG since C > A. Cij{q)is given by 
'At<zi 



aAq) 




B'qt 



A? 

Bi 



A'A^! /q^\ 



qj 



(276) 



above locus FH (i.e., in region 7) and 



a^iq) 




(277) 



below locus FH (i.e., in region 6). 

For small g's (i.e., q±^ <C ^J^, qz ^ £,z^, which corre- 
spo nds to the r egions within the rectangle OJFI in Fig. 
15), < |u((7)p > is given by and C.y (g) is thus 



At iq,T)qi 



{Bq^_+D{T)ql + K{q,T)qiY 

2 9 

mi 



A, {q,T)qlql 



{Bql + D{T)ql+K{lT)qif 



CB^Dq\ 



{Bql + D{T)ql 



K{q,T)qi)^ql 



(278) 



where in the denominator Kq^ is negligible compared to 
Dq^. Now since Cij{q) has three pieces, the crossovers 
are quite complicated. For simplicity, let us divide the 
rectangle OJEI into two subregions separated by locus 
OF, which satisfies Bq^^ ~ Dq\ and hence 



9. = 7^(e^L/e^)^-"-/\i 

SJVL 



(279) 



First we consider the subregion above locus OF (i.e., 
triangle lOF), in which Bq^ dominates Dq\ in the com- 
mon denominator in the Eq. (|278p . In this region, the 
three pieces in the Eq. (|278[) . which correspond to the 



contributions from the the random tilt, the random com- 
pression and the random field respectively, are given re- 
spectively by 



Atqi 




B\t 



'0 

A^ 



qz 

q^_ 

qz 



(280) 
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(281) 



L 



A2 



(e^L/e±) 



9i 



2-r 



(282) 



The crossovers between them is obtained by equating 
them to each other. The function for the crossover (i.e., 
locus DC) between the random compression and random 
field is given by 

:jXo^^NL/t.±) 



TI-Jx^J(^^i/^^)'^'''"'''""''^ (283) 



where we have used the exact scaling relation given in 
the Eq. (|223p : the function for the crossover (i.e., locus 
CE) between the random compression and random tilt is 
given by 



q± 



2 



(284) 



where again we have used the same exact scaling relation 
given in Eq. (|223p : the function for the crossover (i.e., 
locus BC) between the random tilt and random field is 
given by 



These three loci intersect at point C, at which 



(285) 



f- + |i)K-»?t-3 



qz ^ (:,± 



(286) 



It is easy to show that this point is within the triangle 
lOF, as illustrated in Fig. 01 Thus, these three loci divide 
the triangle lOF into three distinct regions with different 
wavevector-dependence of Cij{q) for each region. Specif- 
ically, Cij{q) is dominated by the random compression 
and given by the Eq. (j28ip in region 5, and dominated 
by the random tilt and given by Eq. (|280p in region 4, 



and dominated by the random field and given by the Eq. 
(Pg^ in region 2. 

Now we consider the subregion (i.e., triangle OJF) be- 
low the locus OF, in which Bq^ is dominated by Dq\ in 
the common denominator in the Eq. (|278p . It is easy 
to show that the contribution to Cij{q) from the ran- 
dom compression never dominates in this region. The 
contributions from the random tilt and the random field 
disorder are given by 

'^tqi 



^(r' ^^^^ 



/a) 



2riK-rit 



\2riK-'nt 



(287) 



and 



1 




V Ml 9- 



(288) 



respectively. The crossover (i.e., locus AB) between them 
is 

q^=qfo,a"'~'. (289) 

Cij (q) is dominated by the random field and given by the 
Eq. (|288p in region 1, and dominated by the random tilt 
and given by the Eq. (|287p in region 2. 

Now we discuss the wavevector-dependence of Cij (q) 
in the C-side critical region. For big g's (i.e., gj_ ^ 
or qz » $7^), the wavevector-dependence of Cij (q) is 
the same as that in the A-side critical region. For small 
q's (i.e., qj_ < qz < C!^j{q) is given by 



/ 



a,iq) = 4(g) 



A,, {q'\T)qlqj 



\{7{q',T)ql,+Bql, 
A.' {T)qlql 



K{q',T)ql, 



{j{q'\T)ql,+Bql+K{q\T)qj,y 

{q-',T)ql,ql 
{liq', T)ql, + Bql + K (q', T) q^, 



(290) 

where we have neglected the contribution from the Ax'z' 
disorder (see Hamiltonian (|206p ) since it is not impor- 
tant compared to the contributions from the other disor- 
ders. For convenience, we name the three pieces in the 
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above equation as "A^'-piece", "A^z-piece" and "A^;'- 
piece" , respectively. We find that there are five distinct 
regions in the wavevector space with quahtatively dif- 
ferent wavevector-dependences of Cij (q). Since the az- 
imuthal symmetry about z-axis is now broken due to the 
tilting of the layers, the five regions have to be illustrated 
in a three-dimensional picture, which is shown in Fig. \5\ 
In the following we will show how we derive the five re- 
gions and the wavevector-dependence of Cij (q) for each 
region. 

First we calculate the various crossovers within the 
three planes (i.e., qz'-qx't Qz'-Qs', Qs'-Qx')- Let us start 
with qz'-qx' plane, in which case As/-piece vanishes. One 
crossover in this plane is the scaling locus OG, which di- 
vides the rectangle OHGI into two distinct regions with 
different wavevector-dependences of As', a;' , K and 7. The 
function for OG is given by 



crossover (i.e., locus OA) between A^z-piece and A^/- 
piece are given respectively by 



and 



where 




(294) 



2Cz'/i2Cz'+Vs'). 



(295) 



(296) 



qz'£.z = {qx'^± 



(291) 



The other one is the crossover (i.e., locus OB) between 
A^z-piece and A^j'-piece, which satisfies 



Az'qj, = Ax'ql, 



(292) 



where A^/ is wavevector-indcpendent but temperature- 
dependent and given in (|250p . Ax> has dependences on 
both the temperature and wavevector, which is summa- 
rized in Assuming that OB is in one of the two 
regions separated by OG which is closer to qz'-axis (this 
assumption will be verified as a posterior) , from (j292p we 
obtain 



6 



Vt 



AH ^ ) iizqz') qi. 



A further reorganization leads to 



qz'iz 



where 




'?3/2 



A 



(fe'Ci 



(293) 



A° ^ (^^) A? + A° , 

<^^,z. = 2Cz'/{Kz' +nx')- 

It is easy to verify that the function for the locus OB 
given by the Eq. (|293p is indeed consistent with the as- 
sumption we just made. 

The crossovers in the other two planes can be obtained 
through similar calculations, which will not be repeated. 
We only give the results. In qz'-qs' plane, in which case 
Aa:/-piece vanishes, the scaling locus (i.e., OF) and the 



In qx'-qs' plane, in which case A^'-piece vanishes, coin- 
cidentally, the scaling locus and the crossover between 
Aa;/-piece and As'-piece merge into one (i.e., OD), and 
both are given by the same function: 



iqs'U'^'' 



(297) 



Having found the crossovers in the three planes, now 
we can use them to obtain the crossovers in 3D wavevec- 
tor space. Moving the loci OA and OB along qx' and 
qs' , respectively, we obtain the two boundaries OAC and 
OBC as the trajectories of OA and OB, which meet at 
the locus OC. Thus, OAC and OBC are the crossovers 
between A^z-piece and As/-piece, and A^z-piece and A^;'- 
piece, respectively, and they are also defined by the func- 
tions (|295p and (|293p . respectively, and OC is hence de- 
fined by the the combination of (|295p and (|293p . Like- 
wise, OEF, OEG and ODE are the scaling crossovers, 
which separate distinct regions with different wavevector- 
dependence of Ax' ^s' ^ 7 and K, and arc defined by the 
functions (^Mp . (PST]) and ([Wf)) . respectively. It is easy 
to show that these three crossovers meet at the same 
locus OE, whose function is thus the combination of 
any two of (|^M)l . (POTjl and (^57)) . In addition, ODE 
also serves as the crossover between Aj;/-piece and A^'- 
piece due to the duality of the locus OD. OCE is an- 
other crossover between Aajz-piecc and A^z-piece; how- 
ever, since OCE and ODE are in different regions with 
different wavevector-dependences of Ax' and As' , OCE 
is defined by a function which is different from (|297p : 



{qz'izf 



(298) 



with (ps'x' = [2 - fjK - f}^) /2Cz'- 

After the definitions of the crossovers separating the 
five regions in Fig. [S] become known, the calculation 
of the wavevector-dependence of Cij {q) in each region 
is straightforward. In the region O — ABCH, Cij{q) is 
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dominated by A^/ -piece and given by 



/ 



\ 



\ (7(9', T)ql, + Bql + K (q T) j j 



(299) 



In both of the regions O - AFEC and O - FJDE, C,j {q) 
is dominated by A^/ -piece but with different wavevector- 
dependences, since these two regions are separated by the 
scaUng crossover OEF. In O — AFEC, Cij{q) is given 

by 



A., iq\T)qlql \ 




(300) 



while in 0-FJDE, it is given by 
'A.' {q',T)ql 



K^ql 



^0 9s' 



y \ (<:i- ^2 „6 



^ '^^(e±/e^L)'''-""^'(9s'a)'""-"= 



(301) 



Likewise, in both of the regions O — CEBG and O — 
EDIG, Cij{q) is dominated by A^^'-piece but with dif- 
ferent wavevector-dependences. In — CEGB, it is given 
by 



B2 „4 
,^0 1z 



(302) 



while in O - EDIG, it is 



A(e^j^(^|-) i^^KhLr^'-''^'' 



x{q,,C^)-(^^^+f"<)/<^ 



(303) 



So far we have discussed Cij{q) for small g's (i.e., 
q_L ^ 9z ^ ^7^) in the C-side critical region in 

the transformed reciprocal coordinate system. However, 
what is needed for experiments is the behavior of {q) 
in the lab reciprocal coordinate system. Since the re- 
lation between these two coordinate systems is known, 
the calculation is straightforward but somewhat tedious. 
Instead of listing dj (q) for all qs, which is quite compli- 
cated and far more than necessary, we now discuss Cij (q) 
for special qs, which is simple and adequate for making 
useful experiment predictions. First we consider the case 
in which q is restricted in qz-qx plane. In this case A^/- 
piece vanishes, and Cij (q) reduces to 



Cxxiq) 



A.' (T)g,V 



A,, {q',T)ql,q' 



(7(9 ', T)ql, + Bqi) (7(9 T)ql, + H 



(304) 



The crossover between A^/ and A^,/ has been calculated 
previously and given by Eq. (|293p . which, in terms of q, 
becomes 




A 



{qx - ^qz)^± 



(305) 



This function defines the loci OE and OE' in Fig. [B] 
The coordinates of E and E' are given respectively by 



q^ ^ ±gI^±(6)-^(e±/C^L)^ 



(306) 



and OE and OE' thus connect with EG and E'G' re- 
spectively, which are the crossover between the random 
tilt and random compression fluctuations for big g's (i.e., 
q± ^ S^J^ or qz '3> Cz^)- The crossover between the two 
terms in the common denominator has also been calcu- 
lated and given by Eq. (|29ip , which, in terms of q, is 



qz^z = [(g.+rg,)^^]''^ 



(307) 



38 



It is easy to check that the loci defined by this function 
satisfies {q^l 3> Fq^. Thus this function can be further 
simplified as 

QzL = (q.U)^''^^^' , (308) 

which defines OF and OF' in Fig. [6l These two Loci 
also serve as the scaling crossover separating the re- 
gion for small g's into three distinct ones with differ- 
ent wavevector-dependences of Ax'- Now we calculate 
g-dependences of Cxx{q) in different regions separated 
by these crossovers. In region 4, Cxxift) is dominated by 
A^z-piece and given by 

A.' {T)ql 



qx_ 



in regions 5 and 6, Cxxio) is dominated by A^/-piece and 
given by 

Ax' {q'.T)ql,ql 
mqt 



Bl U^J 4 ^''^'^ 



(^tvl)^ V^jvl/ 1z 



''4 



[Qz^z 



(310) 



in regions 7 and 8, Cxxiq) is also dominated by A^^z-piece 
but given by 

A.' {q\T)ql 



NL 

6 



-(i),,-|-2rj^)/C,, 



(311) 

which is different from Eq. (|310p since regions 5 and 6 
are separated from regions 7 and 8 by the scaling loci. In 
both of the above two equations we approximated g^/ as 
qx, since we have Ig^^l ^ Tqz in regions 5, 6, 7, 8. 

The second special case we consider is when q is along 
qx-Bxis. In this case both A^z-piece and A^j'-piece vanish, 
and Cij {q) is simply given by 

Css{q) 



A? ( U 



\2fiK-flt 



y^NL 

A 



1 



NLl 



2fiK-rit 



(312) 



IX. STABILITY OF TRANSITION AGAINST 
DEFECTS AND ORIENTATIONAL DISORDER 

The stability of the A phase follows from the stability 
of "random field XY model" ; the stability of the C phase 
follows from the stability of "to — 1 smectic". Both 
phases are stable against orientational fluctuations and 
the unbinding of neutral pairs of topological defects (i.e., 
smectic dislocation loops). This implies the existence of 
the phase transition between the two phases. However, 
in order for our theory for the critical behavior to be 
valid, the system also needs to be dislocation-bound and 
orientationally ordered right at the critical point. 

The orientational fluctuations are defined as the mean 
square fluctuations of the layer normal N. Assuming the 
system is orientationally ordered, the dominating orien- 
tational fluctuations can be calculated as 



d'^q 



mil 



i2nr {K{q)qi + BqlY 



(313) 



where Ni_ is the projection of on the _L plane (i.e., x-y 
plane). We have only kept the fluctuations contributed 
by the random tilt disorder, which are the most divergent 
right at the critical point. The (f-dependences of A{q) 
and K{q) are given in Eqs. ((27)) and p2|) . respectively. 
B is g-independent. Then requiring the fluctuations to 
be finite leads to 



2 



< d- 3. 



(314) 



which is one of the conditions for our theory for the crit- 
ical behavior to be valid. 

Now we check if the topological defects are bound right 
at the critical point. The start point of the theory is the 
tilt only Hamiltonian: 



H 



(315) 



where the random field disorder is not included. It can be 
justified that the random field disorder is also irrelevant 
when dislocations are included. Since the model (j315p is 
virtually the same as that for dislocations in smectic A 
phase in an isotropic random environment @ , the theory 
is also the same. Therefore, we will describe the theory 
very briefly. 

When the smectic has a dislocation, the displacement 
field u is no longer single valued. Mathematically this 
can be represented as 



V X VW : 



TO 



with 



rfi{r) 



aM,t{s,)S^ {r - r,{si)) dsi 



(316) 



(317) 
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where Si stands for the whole configuration of the I'th 
dislocation loop, Mi is an integer giving the charge of 
that loop, t{si) is the local tangent of the loop, and fi is 
a parametrization of the path of the loop. Furthermore 
equation (|316p implies 



V • M = 0, 



(318) 



which means that dislocation lines can not end in the 
bulk of the sample. 

To obtain a dislocation Hamiltonian, we need to trace 
over field u which is constrained by equation pi6p . This 
can be done in the following standard way Q . We sepa- 
rate the field v = Vu into 



V = Vd + Sv, 



(319) 



where Vp minimizes Hamiltonian (|315p for a given dislo- 
cation configuration m(r), 5v can be viewed as the fiuctu- 
ation from the ground state. Inserting pi9p into Hamil- 
tonian pi5p . we find that Vp and dv are decoupled due 
to the construction that Vd minimizes Hamiltonian (|315p . 
Thus we obtain the effective model for rn{r). 

Now let us go through the procedure. The Euler- 
Lagrange equation, obtained by minimizing Hamiltonian 
(insD . is 



h = 0. 



Rewriting this in terms of Vd = Vm gives 

1 , 



B 



h = 0. 



(320) 



(321) 



where 
comes 



K/B. In Fourier space, this equation be- 

(322) 



1 



q,v^ + \^q\q±_ ■ vj + —q±_ ■ h{q) = 



and the solution for the constraint (|316p is 
iqxrn 



Vd 



q- 



(323) 



where (p is the smooth elastic distortion around the dis- 
location line. Inserting (|323p into equation (|322p gives 



iqz(l - \^q\)€zijqimj q±-h 



BTn 



(324) 



where we have defined the inverse of the smectic propa- 
gator 



r, = ql + 



(325) 



Inserting (|324p into the solution for Vd and plugging equa- 
tion (|319p into Hamiltonian pisp . we obtain the disloca- 
tion Hamiltonian 



Hd 



d'^q 



Kqj 

2r„ 



mirrij + m ■ a 



(326) 



where = (1 - S,,){1 - *,,)(% - q^ql/ql)-. 
9? + ^^91, and 



qxh {z X q) ■ h 



qz (1 - 



q- 



(327) 



By putting the model on a a simple cubic lattice, the 
partition function can be written as 



Z = 



-S[m] 



where 



m(r) 



(P 

1 

T 



nx{r),ny{r),n,Xr)\ 



Hd[m] 



(328) 

(329) 
(330) 



where the n,;'s are integers, d is the cubic lattice constants 
used in the discretization, and |TO(r)P is the core 

energy term coming from the core of the defect line that 
is not accurately treated by the continuum elastic theory. 

To cope with the constraint V • m = 0, we introduce 
an auxiliary field 0(f), rewriting the partition function 
equation (|328p as 



Z = \{j de{r) J2 e-^["l+''5:.e(F)v.™(F)dVa^ (331) 

r m(r) 

where the sum over m{f) is now unconstrained. 

Then we introduce a dummy vector field A to mediate 
the long-range interaction between defect loops in the 
Hamiltonian (|326p . This is accomplished by rewriting 
the partition function as 

^ = n / d9{r}dA{r) ^ e-^["'^'^l,5(V • A)S{A,) 



(r) 



(332) 



with 



-i V6l(f) + d^[iA{r} + d{r)] 

a 



2T ^ Kq' 



-\A\' 



(333) 



Now the sum over m(r) is recognized to be the "peri- 
odic Gaussian" made by Villain [131 . Then the partition 
function equation (|332p can be rewritten as 



Z 



Y[ J d9{r)dA{f) S{V ■ A)S{A, 



X exp 



E^. 



ad , 



{r + x,)-Oir}-—[A,{r) 



(334) 
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where the 27r-period Villain potential Vp{x) is defined as 

oo 
— 00 

Since Vp{x) has sharper minima for smaller Ec/T, which 
corresponds to higher temperature, raising the tempera- 
ture in the original model is equivalent to lowering tem- 
perature in the dual model equation (|334p . 

Standard universality class arguments imply that the 
model equation (|334p has the same universality class as 
the "soft spin", or Landau-Ginsburg- Wilson, with the 
complex "action" 



E 



ad ,. 1 , « r,-; o-d , . ^, , 

2 [V + — Ma + a)#: • [V - — M„ + a)] 



X4>a + idl^al^ + Ud\(t)a\ 



E 



r, 



(336) 



where ^ is a complex order parameter whose phase is 
0{r). Because of the duality transformation's inversion 
of the temperature axis, the reduced temperature td is a 
monotonically decreasing function of the temperature T 
(of the original dislocation loop model) , which vanishes at 
the mean-field transition temperature of the dislocation 
model p34p . Disorder is included in the model (|336p 
through a(r), which is related to the random tilt field 
h{r) by equation ([327]). 

Because of the duality inversion of the temperature 
axis, the ordered phase of the dual model (|336p corre- 
sponds to the disordered phase (dislocation loops un- 
bound) of the original dislocation model. 

An complete analysis of the dislocation loop unbinding 
transition described by the model p36p is beyond this 
paper. The goal here is to know if a dislocation bound 
state ever exists near the critical region of the smectic A- 
to-C transition. To answer that question, let us check the 
one-loop graphical correction to the dual temperature: 



tR — to 



{d-2)ca^d^ r d'^q ^tqhl 
r2 J {2t:Y q^Tl 



(337) 



It is easy to see that the integral diverges when K and 
At are considered as constants. This implies that Ir 
is always negative and the dual model is always in its 
ordered state, which corresponds to the dislocation un- 
bound state. This implies that the system would have 
become disorder before reaching the critical point, which 
signals either a first order transition, or a reentrant ne- 
matic phase intervening between the A and C phase. 
However, this conclusion only holds within the harmonic 
approximation. In section IIIII we have shown that the 
anharmonic effects are important. A crude way to in- 
clude the effects is to treat K{q) and At(g) anomalous. 
Implementing this in the integral and requiring it to be 
finite, we obtain a restriction on i]x and rjt. 



which is another condition for our theory to be valid. 

We test the two conditions using the values of rjK and 
rit obtained by the e-expansion to O(e^). Unfortunately, 
the two conditions are not satisfied in the physical di- 
mension d = 3, which seems to imply that the second- 
order phase transition is not stable. However, since e = 2 
in d = 3, e-expansion is not quantitatively reliable, and 
therefore whether the second-order phase transition is 
stable in d = 3 remains an open question. 



X. CONCLUSION 

In summary, a theory of smectic A — C phase tran- 
sition in anisotropic disordered media is developed. We 
show that the phase transition can be second-order and 
calculated the critical exponents to the first order in an 
e = 5 — c? expansion. In addition, the elasticity and fluc- 
tuations of this system (at the phase-transition temper- 
ature) are studied. The implications of these results are 
expected to be observable. 



APPENDIX A: 

In this appendix, we review the derivation in reference 
[8| of the numerical estimates of the anomalous elastic 
exponents for the C phase. These are obtained from 
two different e— expansions, based on two different ana- 
lytic continuations of the model to higher dimensions. In 
the "hard" continuation, the higher dimensional (d > 3) 
model is chosen to still have only one "soft" direction for 
all spatial dimensions d, as it does in d = 3. In this case, 
the exponents are given by 



= 2 



2 ' 

Vj = -e + 0(e2), 
Vs' = ^e^ + 0(e~2) 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 



with e = I — d, and obey two exact scaling relations: 



2^7 



7-d. 



-IlK + Vs 



2d- 7. 



(A6) 



2T]t-T]K <0, 



(338) 



In addition, by analytically continuing the problem to 
higher dimensions (d > 3) in a different way (namely, 
by keeping the number of hard directions fixed at 2, as 
opposed to the d = | — e expansion, which is based on 
keeping the number of the soft directions fixed at one), 
this problem was been studied in ^ using a d = 4 — e 
expansion, yielding the following exponents and the exact 
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scaling relation: 





3r 


^0{f), 


(A7) 


V't 






(A8) 


Vs' 






(A9) 


Vs' 


_ Vi 
2 




(AlO) 



4-rf + r)y = 



As argued in [8[, combining the results from the d = 
4 — e and the earlier d = | — e expansion we see that the 
exponents fjfc and ry^' are more consistent between the 
two. This suggests that the most accurate estimate of all 
3 exponents will be obtained by taking rfx and r]s from 
weighted average of the 4 — e and | — e exponents: 



5 

4r?.'(i (4-?) 



(All) 
(A12) 



(the factor of 4 appearing because, a priori, we expect 
the ^ — e expansion to be 4 times as accurate as the 
4 — e expansion, since the errors in both are O(e^), and 
e(= 1) is twice as big as e(= i)), and then obtaining 
rj^ from the exact scaling relation (jA6p in c? = 3 (note 
that (jA10|) reduces to (|A6P in d = 3). One can also 
estimate the errors in the exponents obtained in this way 
as the differences between the weighted averages (IA11[ 
IA12P and the | — e expansion results. Doing all of this 
obtains 



Vs' 



0.50 ±0.03, 
0.26 ±0.12, 
0.132 ± 0.002. 



(A13) 
(A14) 
(A15) 



In addition, inserting 77/4- , 77.^, 773/ into the exact scaling 
relation Eq. gives 



Vx' 



1.37±0.15. 



(A16) 



APPENDIX B: 



During the discussion of the universality class of the 
C phase, we have shown the advantage of writing the 
model in a special coordinate. The relation between this 
coordinate and the lab one is fully determined by the 



parameter F, which is defined by (|205p . Now we show 
that near the critical point, F is also renormalized by 
the critical fluctuations and hence strongly temperature- 
dependent. Using the critical RG, the renormalized F 
can be calculated as 

F = e('^-i)^*F(r), (Bl) 
where the RG is stopped at t* — In(A^j^), and the prefac- 
tor on the right-hand comes from the dimensional rescal- 
ing. To calculate F(^*), it is convenient to reorganize 
(1203) as 



F(r) 



(B2) 



where 33 and are defined by and ([M)) , respectively. 
Previous calculation shows that for large I* (which is true 
near the critical region), gi{t*) — 32e/15 and gsii*) ^ 
e~^^^ , where 773 is defined via 



dgsji) 
di 



(B3) 



According to the linearized flow equation of 93 given by 
(|149p . 773 is given, to the leading order in e = 5 — d, by 



5 

To calculate w'{£*), we rewrite (|197p as 

.93 (r 



gi{£* 



(B4) 



(B5) 



Since gz{i*) / gA{i*) is much less than 1, the above equa- 
tion implies w{l*) / w' {£*) ~ 1. Based on these arguments, 
we obtain from (IB2I1 



F(r 



^ D{t) A ^K{t)' ^ ^ 



where B{i*) and K{£*) can be calculated by integrating 
the RG flow equations (I136i I137p . Plugging the above 
equation into (|Bip . we get 



APPENDIX C: 



TABLE I: Overview of the Elastic Constants and Disorder 



Variances 



Symbols 


Description 


Eq. 


B 


Compression modulus in Hamiltonian Ec[. (|83jl 


B = 


Bo 


K 


Bend modulus in Hamiltonian Eq. (|83p 




At 


Random tilt disorder variance in Hamiltonian Eq. (1831) 


([32|) 


Ac 


Random compression disorder variance in Hamiltonian Eq. (1831) 


1321 


Bo 


Bare value of B 




Ko 


Bare value of K 






Bare value of At 




aL* 


Bare value of Ac 




B 


Compression modulus in Hamiltonian Eq. (12061) 


B = 


Bo 


K 


Bend modulus in Hamiltonian Eq. (12061) 


® 


7 


Tilt modulus along x' in Hamiltonian Eq. (I206II 





A,, 


Random tilt disorder variance along s' in Hamiltonian Eq. (I206II 


m 


A^' 


Random tilt disorder variance along x' in Hamiltonian Eq. (|2U6|I 


® 


A^/ 


Random compression disorder variance in Hamiltonian Eq. (12061) 


A,. = 


= a:, 


Be 


Half-dressed value of B 


Be = 


= Bo 


Kc 


Half-dressed value of K 


(p47l) 


Ic 


Half-dressed value of 7 


(|248l) 




Half-dressed value of A^i 


(|249l) 


A*", 


Half-dressed value of A^i 


(|249l) 




Half-dressed value of A^, 


(|250l) 




TABLE H: Overview of the exponents 








Symbols 


Description 


Eq. 






Anomalous exponent of B 


jys = 






Anomalous exponent of K 


dSS)) 




rjt 


Anomalous exponent of At 


dMI 






Anomalous exponent of Ac 


dSZI 




C 


Anisotropy exponent for the model (I83II 






flK 


Anomalous exponent of K 


m 




fll 


Anomalous exponent of 7 


1121 




fjs' 


Anomalous exponent of A^/ 






fix' 


Anomalous exponent of A^i 






Cx' 


Anisotropy exponent (between g^' and q^i) for the model (12061) 


m 






Anisotropy exponent (between qgi and q^i ) for the model (I206|) 


m 





J 



TABLE HI: Overview of the characteristic lengths 



Symbols 


Description 


Eq. 


A 


Smectic penetration length 


A = ^/A'o/Bo 


SJVi 


Nonlinear crossover lengths 


dMEIl) 




Correlation lengths for the phase transition 


IMD 


a.. 


Crossover lengths between the power-law and broad x-ray scattering patters 






Line widths of the broad x-ray scattering pattern 


(|264M265p 
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